NOV  Tl  IS  SWARNERr  C W CHUANG*  R J HUFF  F30602-75-C-0061 

UNCLASSIFIED  OSURF-710300-1 <ESL> R*DC-TR-7a-»23A ML 


lot* 


UUU  FILE  COPY  AD  AO  62992 


TDMA  TIMING  LOOPS  FOR  HIGH  DATA 
RATE  SYSTEMS 


The  Ohio  State  University  ElectroScience  Laboratory 


W.  G.  Swarner 
C.  W.  Chuang 
R.  J.  Huff 
H.  S.  Eilts 


lApproved  for  public  release;  distribution  uni  I’mitedj 


ROME  AIR  DEVELOPMENT  CENTER 

Air  Force  Systems  Command 

Griffiss  Air  Force  Base,  New  York  13441 


79  01  08  03  3 


- 


This  report  has  been  reviewed  by  the  RADC  Information  Office  (01)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) . At  NTIS  it 
will  be  releasable  to  the  general  public,  including  foreign  nations. 

RADC-TR-78-238  has  been  reviewed  and  is  approved  for  publication. 


APPROVED 


STUART  H.  TALBOT 
Project  Engineer 


APPROVED: 


FRED  I.  DIAMOND,  Technical  Director 
Communications  and  Control  Division 


FOR  THE  COMMANDER: 


JOHN  P.  HUSS 

Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization,  please 
notify  RADC  (DCCR) , Griffiss  AFB  NY  13441.  This  will  assist  us  in  maintaining 
a current  mailing  list. 


Do  not  return  this  copy.  Retain  or  destroy. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 

I REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


REPORT  number 

RABC^-TR-  78-238 

r.  Tl  Til?  (and  Subtitle) 


-n  i- 


|2.  GOVT  ACCESSION  NO.I1.  RECIPIENT’S  CATALOG  NUMBER 


[0  j TDMA  TIMING  LOOPS  FOR  HIGH  DATA  RATE  SYSTEMS  f 


1,  _IYP£  OF  REPORT  A PERIOD  COVERED 


Interim  Report 


6 PERFORMING  ORG  REPORT  NUMBER 

710300-1 (ESL)  

8 CONTRACT  OR  GRANT  NUMBER'S) 


1 — AUTKQRTAi  8 CONTRACT  OR  GRANT  NUMBER'S) 

/PCw.  G.jSwarner.  !,  R.  J.  I Huff  / i 

,'|U  'j  C.  W.jchaung  * ' H.  S.jEilts  ('I  O/F30602-75-C-0061  / L 

— ' ■■  U Q \ ^ ■ ,1  — .. — — 

9 PERFORMING  ORGANIHAJuJl'lfNAME  AND  ADDRESS  '0  P«°<j  u N !V  NUMBE  R$T  ' T *S* 

The  Ohio  State  University  ElectroScience  Laboratory 
Department  of  Electrical  Engineering  ^ 

Columbus  OH  43212  65230907 

11.  CONTROLLING  OFFICE  N AME  AND  ADDRESS  /T~|\  ,aEBOHT"P*Tfi" 

Rome  Air  Development  (DCCR)  ( j | November  1978  , 

Griffiss  AFB  NY  13441  ' -tt~  NUMirETrut-t^TTry^ 


[ II..  -RERORT  OATS  

November  1978  f 


T5  MONITORING  AGENCY  NAME  a ADDRESS  (if  different  from  Controlling  Office)  15.  SECURITY  CLASS  (of  this  report) 


i ! l yjr. 


UNCLASSIFIED 

~is»,  declassifTcation  do*ng r adTn g 

n/aSChedulE 


hi  DISTRIBUTION  STATEMENT  (ot  thi*  Report) 


Approved  for  public  release;  distribution  unlimited 

£ P.  ~J 

17  DISTRIBUTION  STATEMENT  ( o ( the  abstract  entered  in  Block  20,  it  different  Irom  Report) 

s„,  yycur A--". idd HESL ) 


18  SUPPLEMENTARY  NOTES 


RADC  Project  Engineer:  Stuart  H.  Talbot  (DCCR) 


19  KEY  WORDS  (Continue  on  reverse  side  il  necessary  and  identity  by  block  number) 

TDMA  system  timing  Satellite  communications 

Optimum  filters  Synchronization 

Kalman  filters 
Digital  communications 

?0  fContlnue  on  reverse  Side  II  necessary  and  Identify  by  block  number) 

^^Optimum  linear  recursive  filtering  techniques  are  applied  to  achieve  systen 
synchronization  for  a high  data  rate  bit  synchronous  TDMA  satellite  or  relay 
communications  system  employing  two  sampled  data  delay  lock  loops  at  each 
terminal  to  establish  receive  clock  and  transmit  timing  relative  to  a system 
clock  at  the  satellite  or  relay.  It  is  shown  that  a system  employing  an 
optimum  filter  in  the  receive  clock  loop  with  open  loop  coupling  to  a separate 
ranging  loop  for  transmit  timing  can  provide  synchronization  adequate  for  a 
40  Mbit/s  system  containing  highly  maneuverable  airborne  terminals  or  relays. 


DO  1473 


UNCLASSIFIED  / 

SECURITY  CLASSIFICATION  Or  TmiS  P AGE  When  Pei  a Fntere.i 


tlrl  $L>d 


1 


SECURITY  CLASSIFICATION  of  This  PAGE  (When  Data  Entered) 


i 


Filter  models  for  system  timing  are  developed  and  analyzed,  simulation  data  are 
presented  to  verify  and  augment  the  analytical  results,  and  design  curves  and 
equations  applicable  to  receive-only  as  well  as  transmit- receive  synchronization 
for  systems  containing  maneuvering  terminals  or  maneuvering  relays  are  given. 

\ 


PREFACE 


This  report,  OSURF  Report  No. -710300-1  , was  prepared  by  the  Electro- 
Science  Laboratory,  Department  of  Electrical  Engineering,  The  Ohio  State 
University,  Columbus,  Ohio.  The  research  reported  was  funded  by  Rome  Air 
Development  Center,  Griffiss  Air  Force  Base,  New  York,  under  Contract 
F30602-75-C-0061 . Mr.  Stuart  Talbot  is  the  RADC  Program  Monitor. 

The  model  development  and  timing  analyses  of  Section  IV  were 
performed  by  Dr.  Chuang,  the  extended  up-link  analysis  of  Section  V by 
Dr.  Huff,  and  the  simulation  studies  of  Section  VII  by  Mr.  Eilts.  The 
remaining  sections  were  contributed  by  Dr.  Swarner,  who  also  served  as 
task  coordinator. 

The  many  helpful  suggestions  and  guidance  provided  by  Dr.  Huff, 

Mr.  R.  C.  Taylor,  and  others  of  the  OSU  Satellite  Communications  Facility 
relative  to  the  RADC/OSU  TDMA  system  are  acknowledged  and  greatly 
appreciated. 


CONTENTS 


Page 


I INTRODUCTION  1 

II  BACKGROUND  3 

III  THE  BASIC  OPTIMUM  LINEAR  RECURSIVE 

(KALMAN)  FILTER  7 

IV  MODELING  OF  A KALMAN  FILTER  IN  TDM A 

TIMING  SYSTEMS  10 

A.  Introduction  10 

B.  Discrete  Kalman  Filter  for  Clock 

Loop  Timing  12 

C.  Augmented  Discrete  Kalman  Filter  24 

D.  Transient  Response  of  the  Discrete 

Kalman  Filter  35 

E.  Monte  Carlo  Simulation  43 

F.  Up-Link  Synchronization  48 

V UP-LINK  TIMING  - EXTENDED  ANALYSIS  75 

A.  Introduction  and  General  Analysis  75 

B.  Stationary-Satellite  Maneuvering-Terminal  Case  91 

C.  Stationary-Terminal  Maneuvering-Satel 1 i te  Case  97 

VI  APPLICATIONS  TO  SYSTEM  DESIGN  111 

A.  Clock  Loop  Timing  111 

B.  Up-Link  Timing  130 

C.  Application  of  the  Extended  Up-Link 

Timing  Analysis  140 

VII  SIMULATION  STUDIES  174 

A.  Transient  Response  174 

B.  Steady  State  Zero  Motion  Response  179 

VIII  CONCLUSIONS  184 

REFERENCES  186 


SECTION  I 
INTRODUCTION 


This  report  describes  the  application  of  linear  optimum  recursive 
filtering  to  upgrade  system  timing  for  a bit-synchronous  time  division 
multiple  access  (TDMA)  relay  communications  system.  This  system  was 
previously  developed  at  The  Ohio  State  University  (OSU)  ElectroScience 
Laboratory  [1-6]  under  contracts  v/ith  the  Rome  Air  Development  Center 
(RADC)  of  the  United  States  Air  Force.  Greatly  increased  signaling 
rates  and  improved  tracking  of  highly  maneuverable  airborne  terminals  or 
relays  are  thus  achievable  while  retaining  the  advantages  of  efficiency, 
flexibility,  and  anti-jam  protection  capabilities  inherent  in  this  system. 

In  Section  II  background  information  relative  to  the  RADC/OSU  TDMA 
system  is  given,  and  in  Section  III  a brief  summary  of  the  pertinent 
results  of  basic  Kalman  filter  theory  [7,8]  are  presented  as  a preface 
to  the  analyses  of  Sections  IV  and  V. 

In  Section  IV,  a two-dimensional  model  for  the  clock  loop  timing 
filter  is  first  developed  and  transient  and  steady  state  solutions  are 
obtained  following  essentially  the  method  employed  by  Friedland  [9].  A 
three-dimensional  augmented  model  is  then  developed  employing  an 
acceleration  correlation  coefficient  similar  to  that  used  by  Singer  and 
Behnke  [10]  for  whitening  the  maneuver  noise.  Closed-form  steady  state, 
as  wefl  as  transient  solutions,  for  this  model  are  obtained.  Transient 
response  is  then  studied  utilizing  computer  iteration  techniques,  and 
results  of  Monte  Carlo  simulation  studies  to  check  the  validity  of  the 
models  and  analytical  results  are  presented.  Finally,  the  problem  of 
up-link  synchronization  is  studied  using  a method  of  open  loop  cor- 
rections from  the  clock  loop  to  the  ranging  loop  to  augment  the  closed  loop 
transmit  timing  corrections.  Expressions  for  the  clock  loop  and  ranging 
loop  timing  errors  are  derived  and  closed-form  steady-state  solutions 
for  the  transmit  timing  error  variances  for  a system  containing  a maneuver- 
ing relay  and  for  a system  containing  maneuvering  terminals  are  obtained 
and  evaluated. 

Section  V presents  a more  complete  "second  generation"  analysis  of 
up-link  timing  utilizing  a modified  open  loop  feedback  scheme  to  achieve 
somewhat  lower  error  variance  under  receiver  noise  limited  (as  opposed  to 
maneuver  limited)  conditions.  This  analysis  also  features  retention  of 
additional  cross  correlation  terms,  neglected  in  the  previous  analysis, 
which  should  lead  to  more  accurate  analytical  results  at  a cost  of  increased 
complexity  of  the  resulting  expressions  and  a corresponding  increase  in 
computational  effort  required  to  obtain  numerical  data. 

Applications  of  these  results  to  system  design  are  considered  in 
Section  VI.  The  form  of  the  steady  state  solutions  of  Sections  IV  and 
V,  although  mathematically  concise,  are  not  particularly  amenable  to 
direct  application  in  system  design;  hence  alternative  forms  are 
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developed  and  a series  of  system  design  curves  giving  system  parameters 
as  a function  of  known,  or  specified,  signal  and  vehicle  maneuverability 
statistics  and  required  timing  accuracy  are  presented.  Comparisons 
between  the  approximate  solutions  for  up-link  timing  of  Section  IV  and 
the  more  exact  solutions  of  Section  V are  also  given. 

In  Section  VII,  typical  examples  of  computer  simulation  results  are 
given  showing  system  transient  response  and  steady  state  performance 
for  different  signal  and  noise  parameters,  and  a comparison  of  steady 
state  error  variance  with  analytical  results  is  given.  Conclusions 
and  recommendations  are  summarized  in  Section  VIII. 


SECTION  II 
BACKGROUND 


Techniques  for  implementing  time  division  multiple  access  (TDMA) 
satellite  communications  systems,  particularly  those  containing  a 
large  number  of  small  non-stationary  terminals,  e.g.,  aircraft,  have 
previously  been  developed  and  successfully  demonstrated,  both  theoreti- 
cally and  experimentally,  using  prototype  modems  and  a satellite 
simulator  with  adaptive  spatial  processor  developed  and  constructed  at 
this  laboratory  [1-6]. 

The  signaling  format  for  this  system  is  shown  in  Figure  1 [5], 

The  time  continuum  is  divided  into  non-overlapping  intervals  or  slots, 
each  of  which  is  (normally)  allocated  for  the  relaying  of  signal  from 
no  more  than  one  terminal  at  a time.  The  slots  are  defined  with  respect 
to  the  time  base  of  the  netv/ork  clock  signal  (NCS)  present  on  the 
satellite  down-link  once  each  subframe.  At  each  user  terminal  the 
time-base  of  a locally-generated  signal  (clock)  is  aligned  with  the 
time  base  of  the  received  NCS  to  establish  a local  receive  clock.  In 
turn,  a transmit  clock  is  timed  so  that  pulses  transmitted  by  the  terminal 
occupy  assigned  time  slots  on  arriving  at  the  satellite.  Proper  trans- 
mit clock  timing  is  obtained,  during  the  link-range  slots,  by  estimating 
the  error  in  arrival  time  of  pulses  transmitted  by  the  terminal  as  they 
are  received  on  the  down-link  relative  to  the  local  receive  clock. 

The  principal  features  of  this  system  are: 

(1)  Only  one  signal  is  incident  on  the  satellite  at  a time  so 
that  small  terminals  do  not  have  to  compete  for  down-link  power 
with  much  larger  ground  based  terminals  which  may  also  be  using 

the  system. 

(2)  A flexible  demand  assignment  signaling  format  permits 
efficient  use  of  channel  capacity  for  a large  number  of  both  small 
and  large  terminals  at  various  data  rates  compatible  with  their 
needs  and  capabilities. 

(3)  Differential  detection  is  used  (with  2<}>  or  4$  DPSK 
modulation)  to  minimize  signal  coherence  requirements  so  that  an 
adaptive  spatial  processor  (adaptive  array  antenna)  can  be  ef- 
fectively used  at  the  satellite  for  up-link  protection  from 
interference  and  jamming. 

(4)  Bit  synchronous  timing  of  both  up-link  and  down-link 
signals  relative  to  a system  clock  at  the  satellite,  and  the  use 
of  differential  detection,  essentially  eliminate  the  need  for 
preambles  and  guard  space  between  data  bursts  thus  permitting  ef- 
ficient signal  switching  and  adaptive  array  operation  for  small  as 
well  as  large  data  bursts. 
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Figure  1.  The  TDMA  modem  signaling  formats. 


A moderate  amount  of  spectrum  spreading,  via  a suitable  PN  code, 
is  used  as  a basis  for  identification  of  desired  signals  for  spatial 
processing  and  also  provides  additional  processing  gain  for  inter- 
ference reduction. 

To  obtain  reliable  operation  the  timing  of  all  received  signals 
relative  to  the  receive  clock  must  be  accurately  maintained  at  each 
terminal.  Practical  timing  error  limits  are  about  5 percent  of  the 
spectrum  spreading  code  chip  width  to  avoid  a significant  increase  in 
error  rate  [1].  In  the  present  system  this  requires  accurate  control 
of  the  transmit  timing  as  well  as  the  receive  clock  timing  relative  to 
the  system  clock  at  the  satellite  and  is  accomplished  by  means  of  two 
sampled  data  delay  lock  loops  at  each  terminal.  The  first  loop  phase 
locks  the  local  receive  clock  to  timing  pulses  received  on  the  down-link 
from  the  satellite  once  each  subframe.  These  timing  pulses  may  either 
originate  at  the  satellite  - which  is  especially  advantageous  when  an 
adaptive  array  is  employed  at  the  satellite,  or  may  be  relayed  by  the 
satellite  from  a designated  network  control  terminal.  With  the 
receive  clock  timing  thus  established  the  terminal  is  able  to  receive 
all  properly-timed  transmitted  signals  relayed  by  the  satellite.  The 
second  loop  then  adjusts  the  terminal  transmit  clock  so  that  a ranging 
pulse  transmitted  by  the  terminal  once  each  frame  is  received  on  the 
down-link  properly  synchronized  with  the  receive  clock  timing  previously 
established.  These  delay  lock  loops  operate  linearly  provided  that  the 
timing  error  does  not  exceed  ±1/2  chip.  Separate  coarse  ranging  (search) 
circuits  (not  considered  in  this  report)  are  used  for  initial  acquisition, 
or  after  loss  of  lock,  to  bring  the  error  within  this  operating  range. 

The  prototype  system,  designed  to  demonstrate  concepts  and 
feasibility  and  to  obtain  experimental  data  compatible  with  the  small 
terminal  concept,  employed  a relatively  low  code  chip  rate  of  175.2  Kbps 
for  the  lower  rate  format,  and  1.4016  Mbps  for  the  higher  rate  format. 

The  required  timing  accuracy  for  this  system,  for  both  transmit  and 
receive  clock  timing,  was  readily  achieved  using  simple  averaging 
filters  (no  memory)  in  the  delay  lock  loops.  As  the  data  rate  and/or 
the  terminal  maneuverability  is  increased,  however,  absolute  timing  accuracy 
must  also  increase  in  order  to  retain  the  same  relative  accuracy  and 
corresponding  error  rate.  The  techniques  developed  in  this  report  should 
provide  timing  accuracy  sufficient  for  similar  systems  operating  at  4G 
Mbit/s,  or  greater,  with  highly  maneuverable  terminals  or  relays  having 
accelerations  on  the  order  of  5 g. 

The  method  of  improving  timing  accuracy  developed  in  this  report 
employs  an  optimum  linear  recursive  (Kalman)  filter  in  the  terminal 
receive  clock  loop  to  derive  an  optimum  estimate  of  the  system  timing, 

T,  (time  delay  proportional  to  satellite  range)  from  the  noise-corrupted 
clock  pulses  received  from  the  satellite  by  the  maneuvering  terminal. 

Optimum  estimates  of  delay  rate  t (velocity)  and  T (acceleration)  are 
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also  derived  by  this  filter  and  are  fed  as  open  loop  corrections  to 
the  transmit  timing  loop  in  a manner  similar  to  that  used,  for  T only, 
in  the  original  system  [1].  Thus  the  prediction  capability  of  the 
clock  loop  filter  is  utilized  in  the  ranging  (transmit  timing)  loop  as 
well.  (Unfortunately,  timing  errors  due  to  clock  oscillator  drift  are 
essentially  doubled  when  the  open  loop  correction  scheme  is  used.) 

Note  that  to  prevent  instability,  the  rate  at  which  new  range 
estimates  can  be  derived  by  the  closed  loop  ranging  loop  is  limited 
by  the  round  trip  path  delay  between  the  terminal  and  the  relay 
(approximately  0.25  sec.  for  a synchronous  satellite  relay),  while  no 
such  restriction  applies  to  the  system  clock  pulse  rate  and  terminal 
clock  loop.  Hence  for  a satellite  communications  system  these  open  loop 
corrections  can  be  made  at  a much  higher  rate  (several  times  per 
subframe)  than  can  the  closed  loop  range  estimates  derived  by  the  ranging 
loop  (once  per  frame). 

In  effect  then,  the  clock  loop  provides  the  tracking  capability 
for  transmit  timing  as  well  as  receive  timing,  except  for  a constant 
range  offset.  Only  this  constant  offset,  plus  corrections  for  in- 
accuracies in  the  open  loop  inputs  and  for  oscillator  drift  needs  to 
be  determined  by  the  ranging  loop,  and  this  can  be  done  effectively  at 
a slower  rate  with  a less  sophisticated  circuit.  For  this  reason,  in 
the  present  analysis,  the  original  averaging  filter  is  retained  in 
the  ranging  loop.  However,  some  improvement  in  transmit  timing  accuracy 
should  be  achievable  by  applying  an  optimum  filter  in  this  loop  also  to 
obtain  a better  estimate  of  the  offset  and  drift  parameters  and  this  may 
be  the  subject  of  future  studies. 
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SECTION  III 

THE  BASIC  OPTIMUM  LINEAR  RECURSIVE  ( KALMAN) FILTER 


The  basic  theory  of  the  discrete  optimum  linear  recursive  filter, 
often  called  the  Kalman  filter,  is  well  known  and  adequately  documented 
elsewhere  (see,  for  example.  Sage  and  Melsa  [7]  and  Nahi  [8])  and 
will  not  be  repeated  here.  However,  some  pertinent  results  and  im- 
plications of  the  basic  analysis  which  are  applicable  to  the  model 
development  and  subsequent  analytical  solutions  of  Sections  IV  and  V 
are  summarized  here  for  convenience. 

A convenient  form  of  the  basic  Kalman  filter  for  application  to 
the  present  problem  is  given  in  Figure  2.  The  purpose  of  this  filter 
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Figure  2.  Basic  Kalman  filter. 


is  to  derive  an  optimum  filtered  estimate  Zn  of  the  system  state  vector 
Zp  following  arrival  of  each  noise-corrupted  measured  data  point  Xn. 

As  shown  in  Figure  2,  this  is  accomplished  by  summing  two  components: 
(1)  the  predicted  estimate  Zn,  and  (2)  the  innovation  c^n  multiplied  by 
the  Kalman  gain  vector  Kn. 

The  predicted  estimate  ^n,  which  provides  the  memory  of  "learned" 
component^of  each  new  estimate,  is  derived  from  the  previous  filtered 
estimate  Zn_i  via  the  transition  matrix  41.  Thus  the  predicted  estimate 
is  up-dated  and  recycled  with  the  receipt  of  every  data  point  so  that, 
at  any  given  time,  this  estimate  contains  information,  with  optimum 
weighting,  on  all  previous  data  points  collected  up  to  that  time. 
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The  innovation^  contains  the  new  information  supplied  by  the 
current  measured  data  point  and  is  obtained  by  subtracting,  from  this 
measured  data,  the  predicted  measurement  which  is  the  best  estimate 
of  what  the  measurement  should  have  been,  based  on  all  previous  data 
collected  to  date.  The  predicted  measurement  is  obtained  from  the 
predicted  estimate  via  the  H-matrix  which  allows  for  the  fact  that 
the  measurement  variables  Xn  may  differ  from  the  system  state  variables 
Zn.  The  innovation  therefore  contains  only  the  random  part  (white  noise) 
of  the  measured  data  and  consists  of  (for  the  present  example)  the 
random  component  of  the  effects  of  terminal  maneuvers  and  the  measure- 
ment noise. 

The  Kalman  gain  vector  provides  the  optimum  weighting  between  the 
innovation  and  the  predicted  estimate  and  normally  varies  (decreases) 
with  each  subsequent  data  point  after  initial  turn-on  or  acquisition 
until  steady  state  is  reached  asymptotically,  at  which  time  the 
adaptive  Kalman  filter  is  equivalent  to  the  steady  state  Wiener  filter. 

The  Kalman  gain  vector  is  computed  for  each  specific  model  (see  Section 
IV),  and  its  value,  for  each  iteration  n,  depends  upon  the  prior 
statistics  assumed;  i.e.,  the  measurement  noise  arid  terminal  maneuver 
noise  statistics  which  must  be  supplied  externally.  Note  that  the  gain  in 
no  way  depends  upon  the  measured  data  (a  common  misconception  per- 
petrated by  the  often  used  term  - adaptive  filter)  and  hence  can  be 
pre-computed  and  stored  if  desired,  rather  than  being  computed  iteratively 
concurrent  with  the  measurements.  However,  if,  in  addition  to  the  basic 
Kalman  filter,  a means  is  provided  for  measuring  and  upgrading  the 
prior  statistics  as  a function  of  time,  these  may  be  used  in  the 
iterative  gain  calculations  to  provide  true  adaptability. 

The  error  variances,  i.e.,  the  variance  of  the  difference  between 
the  actual  value  of  each  state  variable  and  its  estimated  value,  are 
also  normally  calculated  (since  these  are  used  in  calculating  the  gains) 
and  provide  a measure  of  filter  performance.  These  calculated  variances, 
however,  like  the  gain  calculations,  are  accurate  and  optimum  only  to 
the  extent  that  the  actual  measurement  noise  and  terminal  maneuver 
statistics  agree  with  the  prior  statistics  assumed  for  the  calculation 
and,  of  course,  to  the  degree  that  the  model  accurately  describes  the 
actual  physical  process.  Hence  computer  simulation  or  actual  measure- 
ments are  highly  desirable  in  addition  to  the  filter  analysis  to  validate 
the  results  obtained. 

As  is  evident  from  Figure  2,  increasing  the  Kalman  gain  vector 
increases  the  relative  input  from  the  current  measurement  and  hence 
increases  the  capability  of  the  filter  to  track  a maneuvering  terminal. 

It  also  increases  the  contribution  from  measurement  noise,  hov/ever, 
and  so  it  is  not  surprising  that  the  optimum  gain  setting  is  found  to 
depend  upon  the  ratio  oa/of,  where  aa  and  a ^ are  the  standard  deviations 
of  the  maneuver  "noise"  and  measurement  noise,  respectively.  Increasing 
this  ratio  increases  the  optimum  gain  setting,  and  conversely.  In 
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particular,  if  we  let  the  maneuver  statistics  (oa)  approach  zero  (i.e., 
stationary  terminals  and  relay)  then  the  gain  vector  asymptotically 
approaches  zero  with  increasing  time  (since  the  state  vector  becomes 
completely  predictable,  and  hence  the  gain  can  be  reduced  to  zero 
to  reduce  the  input  noise),  and  the  filter  becomes  decoupled  from  the 
measurement.  This  situation  must  obviously  be  avoided  since,  if  this 
happens,  even  the  smallest  error  in  the  modeling  and/or  measurement 
process  will  lead  to  very  large  errors  in  the  predictions  even  though 
the  predicted  error  variance  will  approach  zero  along  with  the  gain,  a 
condition  known  as  divergence.  This  problem  is  easily  avoided  by 
assuming  a sufficiently  large  value  of  oa  in  the  modelling  process. 

Then  if  the  actual  maneuver  statistics  fall  below  this  value  the 
filter  performance,  although  not  truly  optimum,  will  generally  be 
adequate,  and  divergence  will  not  occur. 

The  principal  equations  associated  with  the  basic  Kalman  filter  of 
Figure  2 are  summarized  below  for  convenience. 


A, 

Zn  = 


(filtered  estimate) 

(1) 

n-1 

(predicted  estimate) 

(2) 

% 

(3) 

- H Zn 

(innovation) 

where  Zp  is  the  system  state  vector  modeled  recursively  as 


-n+1 


z + r u 
n 'tr 


(4) 


and  the  measurement  is  given  by 


X = H Z + 5 . 
i n 'vn 


(5) 


In  the  above  equations,  and  are  white  maneuver  noise  and 
measurement  noise,  respecti vely.  The  state  vector  Zn,  state  transition 
matrix  i,  and  the  matrices  H and  r are  model  dependent  and  are  developed 
for  two  specific  models  in  the  following  section,  where  expressions  for 


the  Kalman  gain  vector 
are  also  derived. 


Kn  and  the  associated  error  covariance  matrices 


SECTION  IV 

MODELING  OF  A KALMAN  FILTER  IN  TDMA  TIMING  SYSTEMS 


A.  Introduction 


The  terminal  receive  clock  timing  system  of  the  TDMA  modems 
developed  at  The  Ohio  State  University  ElectroScience  Laboratory 
employs  a sampled  data  delay  lock  loop  (SDDLL)  as  shown  in  Figure  3a 
[1].  When  in  lock,  the  sampled  error  voltage  versus  timing  error 
characteristic  (Es  vs  e)  is  linear  for  |e]  < a/2  where  A is  the  chip 
length  of  the  PN  code.  This  sets  an  upper  limit  on  the  timing  error 
for  linear  operation  of  the  SDDLL  at  half  the  chip  length.  The  discrete 
loop  filter  in  the  SDDLL  is  a simple  averaging  filter  in  which  N=4 
samples  are  currently  averaged  to  obtain  an  output  error  voltage  sample 
Es  which  then  is  used  to  adjust  the  clock  rate  to  reduce  the  error.  A 
simplified  model  of  the  SDDLL  including  the  averaging  filter  is  given 
in  Figure  3b. 

Values  of  the  chip  length  A and  the  timing  correction  period  Tf 
(4/subframe  rate)  for  the  current  system  in  which  4 pulses  are  averaged 
are  5.7  us  and  107  ms,  respecti vely,  for  the  lower  rate  format;  and 
0.71  usee  and  13.3  ms,  respectively,  for  the  higher  rate  format  giving 
a ratio  Tf/A  = 1.9  • 10^  for  both  data  formats.  If  we  allow  a timing 
error  of  .05a  for  acceptable  performance  (to  avoid  a significant  increase 
in  error  rate;  then  a range  rate  up  to  2.6  usec/sec  (%18G0  mi/hr  relative 
terminal  velocity-assuming  no  oscillator  drift  or  noise  contributions) 
could  be  accommodated.  In  the  ranging  (transmit  timing)  loop,  however, 
where  the  minimum  correction  time  is  limited  by  the  round  trip  path  delay, 
the  correction  time  used  is  Tp  = frame  rate  = 0.427  sec.  yielding,  for 
the  same  error  criteria  as  above,  a maximum  range  rate  of  0.67  usec/sec 
(450  mi/hr)  for  the  lower  rate  format  and  0.084  usec/sec  (56  mi/hr.) 
for  the  higher  rate  format.  Thus  the  need  for  open  loop  corrections 
to  the  ranging  loop  from  the  receive  clock  loop  (currently  provided  by  the 
cross  strap  - see  Section  II)  is  clearly  indicated  when  maneuvering 
terminals  and  high  code  chip  rates  are  involved. 

We  now  desire  to  increase  the  code  chip  rate  to  approximately 
40  Mbps,  or  greater,  from  the  current  values  of  175  Kbps  and  1.4  Mbps 
for  the  lower  and  higher  rate  formats,  respectively,  and  to  allow  for 
tracking  of  rapidly  maneuvering  terminals.  For  a relative  timing  error 
of  0.05a  an  absolute  timing  accuracy  on  the  order  of  1 nsec  will  now 
be  required. 

In  the  absence  of  measurement  noise,  the  required  tracking  accuracy 
could  be  achieved  simply  by  scaling  the  sampling  interval  by  the  same 
ratio  as  the  chip  length,  i.e.,  retaining  the  same  TDMA  format  as 
before,  as  measured  in  chips,  since  the  relative  timing  error  due  to 
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(b)  MATHEMATICAL  MODEL 

Figure  3.  The  sampled  data  delay-lock  loop  (SDDLL) . 


Ol-«> 


vehicle  maneuvers  depends  on  the  ratio  Tf/A.  The  relative  error  due  to 
measurement  noise,  however,  is  (for  E/N0  > 10,  UQ  = constant)  approxi- 
mately proportional  to  (1/E)'/^  (see  Equation  (48)  of  Reference  1), 
where  E is  the  timing  pulse  energy  and  N0  is  the  measurement  noise 
density;  and  E can  be  held  constant,  with  decreased  a,  only  by  increasing 
the  received  power  pr  and/or  increasing  the  number  of  chips  M per 
timing  pulse  (i.e.,  using  a less  efficient  TOMA  format).  Thus  the  need 
for  a more  efficient  timing  system  for  higher  data  rate  TDMA  systems 
of  this  type  is  evident.  As  mentioned  earlier  in  Section  II  this  is 
accomplished  by  incorporating  an  optimum  linear  (Kalman)  filter  into 
the  receive  clock  timing  loop  and  using  the  prediction  and  tracking 
capabilities  of  this  filter  also,  by  open  loop  coupling,  in  the 
ranging  (transmit  timing)  loop.  Appropriate  models  for  this  filter 
wi 1 1 now  be  derived. 


B.  Discrete  Kalman  Filter  for  Clock  Loop  Timing 

A two-dimensional  model  of  the  Kalman  filter  applicable  to  the 
clock  loop  timing  SDDLL  is  developed  in  this  section. 

Let  T be  the  delay  of  the  timing  signal  (proportional  to  the  range 
between  the  terminal  and  the  satellite),  T be  the  delay  rate  and  T 
the  rate  of  change  of  the  delay  rate,  and  let  the  observation  interval 
be  Tf.  The  observation  is  corrupted  by  a Gaussian  white  noise  i.e.. 


X(n)  = T(n)  + C(n) 


(6) 


where  n represents  the  nth  sample,  and  the  expected  values  of  the  measure- 


ment  noise  are  given  by 

EU(n)]  = 0 

(7) 

E[^(n)]  = = const. 

(8) 

E[c(n) s(k) ] = 0 (n/k)  . 

(9) 

The  relative  motion  between  the  terminal 
proximately  described  by 

and  the  satellite  can  be  ap- 

T(n+1 ) = T(n)  + T(n)Tf  + \ T(n)T^ 

(10) 

T(n+1)  = T(n)  + T(n)  T 

(ID 
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If  T(n)  can  be  modeled  by  white  noise,  then 


E[T(n) ] = 0 

02) 

E[T^(n)]  = = const. 

a 

03) 

E[T(n)T(k)]  = 0 (n/k)  . 

(14) 

Writing  (10),  (11)  and  (6)  in  matrix  form,  we  have 

Z(n+1 ) = * Z(n)  + rf(n) 

05) 

where 


Z(n) 


r 


and 


T(n) 

T(n) 


0 1 


X ( n ) = H Z(n)  + c(n) 


where 


(16) 

(17) 

(18) 

(19) 


H = [1  0].  (20) 

For  the  specific  appl i catjon..of  the  Kalman  filter  in  the  SDDLL,  it 
is  more  convenient  to  write  T,  T,  T,  and  oa  in  units  of  the  chip 
length  a.  The  observation  interval  Tf , which  is  controlled  by  the 
local  clock,  is  not  constant  because  of  the  non-zero  delay  rate  f. 

However,  as  shown  below,  the  observation  interval  Tf  in  (10)  and  (11 ) 
can  be  replaced  by  its  constant  nominal  value  with  negligible  error. 

Assume  that  the  code  frequency  of  the  timing  signal  is  40. MHz.  Then 
the  chip  length  a is  equal  to  25  nsec.  Also  assume  that  T = 3 usec/sec  = 
120  a/s ec  (about  three  times  the  speed  of  sound)  and  the  nominal  ob- 
servation interval  = .01  sec.  The  change  in  the  observation  interval  from 
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one  observation  to  the  next  is  only  30  nsec.  If  this  change  in  the 
observation  interval  is  included  in  (10)  and  (11),  it  v/i  1 1 only  con- 
tribute t x 30  nsec  = 3.6  x 10““  A,  which  is  negligible.  Therefore,  the 
observation  interval  Tf  will  be  considered  constant. 

Following  standard  procedures  for  the  Kalman  filter  (see  Chapter 
III)  we  obtain  the  optimum  estimate  of  Z(n), 


Z(n)  = ^(n)  + K(n)  [X(n)  - H £(n)] 

(21) 

^(n)  = «>  Z ( n- 1 ) 

(22) 

where 


Z(n)  = 


is  the  optimum  estimate  of  Z(n)  after  the  measurement  X(n)  is  processed 
(i.e. , the  filtered  estimate),  and 


T(n) 
t ( n ) 


(23) 


2(n)  = 


l'(n) 

_t(n) 


(24) 


is  the  optimum  estimate  of  Z(n)  before  the  measurement  X(n)  is  processed 
(the  predicted  estimate).  The  gain  matrix  K(n)  is  given  by 

K(n)  = P(n)  HT[H  P(n)  HT  + R]'1  (25) 


where 

R = a*  (26) 

and 

ftn)  = E[(^(n)  - Z(n))(2(n)  - Z(n))T]  (27) 


is  the  covariance  matrix  of  estimation  error  prior  to  processing  X(n). 
The  covariance  matrix  P(n)  is  computed  recursively  using  the  variance 
equation 


^(n+l)  = *P(n)$T  +rqrT 


(28) 


wi  th 


P(n)  = [ I-K( n)  H]  P(n)  (29) 

where 

P(n)  = E[( Z(n)  - Z(n))(Z(n)  - Z(n))T]  (30) 

is  the  covariance  matrix  of  estimation  error  after  processing  X(n), 

o ■ i on 

is  the  variance  of  the  random  acceleration,  and  superscript  T indicates 
the  transpose  of  the  matrix.  Using  the  matrix  inversion  lemma  [7],  we 
obtain  from  (29), 

P_1(n)  = p’1  (n)  + hY1  H.  (32) 

Multiplying  (32)  by  P(n)l  )P(n),  we  obtain 


?(n)  = .P(n)  + P(n)HTR_1  H P(n). 


(33) 


By  comparing  the  above  equation  with  (29),  we  have 
K( n)  = P(n)  HT  R'1  . 


(34) 


The  initial  estimation  of  Z (after  two  data 
is  given  by 


Z(2)  = 


X(2) 

f-  (X(2)-X(l)) 
'f 


points  have  been  received) 


(35) 


where  the  values  of  X(l)  and  X(2)  are  obtained  implicitly  (outside  the 
SUDLL)  from  the  1st  and  2nd  range  delay  measurements,  respectively.  Then 
the  initial  value  of  the  covariance  matrix  of  estimation  error  is 


P(2) 


aVTf 


la2  T2 


2 ,T2 


f + 2oVTf 


(36) 
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The  transient  behavior  is  best  studied  using  numerical  methods 
and  will  be  given  in  section  D.  The  steady-state  behavior  is  given  in 
the  following. 


In  the  steady  state,  (28)  and  (29),  using  (25)  reduce  to 


and 


^ - T T 

p = * p + rqr1 


" 'X.'VT  % T -1  % 

P = P - PH  [HPH  +R]  1 H P 


(37) 

(38) 


Combining  (37)  and  (38), 

(p-rqrT)($_1 )T  = p-pht[hpht+r]'1hp. 


Let 


% 

P 


% 

P2 

% 

P, 


(39) 


(40) 


then 


HPHT+R=[1  0] 


r 

P1  P2 

1 

% 'V 

h P_3_ 

_0_ 

2 -v  ? 
+ a;  = P,  + at 


(41) 


^ht[h^ht+r]hp  = 


1 


2 

P1+0C 


P1 

p 

K2 

r 

[1  0] 

'V 

P1 

f 

*2 

1 

p 

P f* 

yr  2 

P2 

p2 

_o_ 

'V, 

[_P2 

PL 

% 2 
P1+°C 

L!lP2 

p2 

2_ 

and 


^ ^ Tr-  'V  T -.-1  'V  1 

P-PH'[HPH'+R]  1 HP  = 


P1+0C 


VI 

-oe 

ro 

Q 

hi 

pip: 

\o  ^ ? 
2 T? 


(42) 


(43) 


We  also  have 
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Define 


% 

A = !i 

A 2 


% 

Po 


B = 


°£°aTf 


and 


r = 


4o? 

°7? 


Then  from  (50): 

A = B2  - 1 . 


(48)  x if  + (49)  x -£  : 


A-iB  + 4 = A + B 4 

r T 1 +A  T+A  r 


From  (55): 


A2+4(l+A)/r2 

2+A 


Substituting  (56)  into  (54), 


.4,4  „ 2,3 


2\  ,2 


r A -8r  A +(16-72r  )A  +(32-128r  )A+(16-64r  ) = 0 


or 


[r2A2-4(l+2r)A-4(l+2r)][r2A2-4(l-2r)A-4(l-2r)]  = 0. 


Therefore 


, _ + */l'+2r  (/l+2r  i I)2 

M 2 

r 


L 
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(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 


P 


or 


« _ + /U2r  (/U2r  ± ^)2 

A ^ * (60) 

The  P?r^ive  sign  °Ptions  must  be  chosen  since  A > 0 and,  the  second  solu- 
tion (60)  must  be  discarded  because  it  makes  P-|P 3 - p|  < o.  Therefore: 

A = v 1 +2 r + 1 )2 


Final ly  we  obtain 


J_  _ *l+2r  (/l+2r  + ^)c 


a a T£ 
5 a f 


(/l+2r  + 1)‘ 
2r 


£ - 

-fj  = \ (^  + i) 

°aTf 


where 


The  ratio  r is  a dimensionless  parameter  which  could  be  regarded  as  a 
form  of  noise  to  signal  ratio. 

From  (38) 


^ P p p 

1 P2  P1  p2 

D /v 

P2  P3  + T P2  P 3 


where 


_ a-  % 

0 ■ P1P3-PI 


Therefore 
P 


-i  = 'r^r'  ( ^ ^ +2r  - l)2 


a a T, 

5 a f 


_ (/l +2r  - 1)‘ 

2r 


h 7 = ? (/J7^  - ]) 


a T, 
a f 


Then  from  (34) 

/ 

'T 


K 


LKf 


= P hV1 


V»c 


UV0?. 


(67) 


(68) 


(69) 


(70) 


(71) 


Therefore  the  steady  state  gain  factors  are  given  by 


and 


M+2r  (i 


l+2r  - 
T~ 


1 )‘ 


(72) 


= 2(rT+2r  - I)2 


V 


(73) 


/v  2 ^ 2 

Graphs  of  P|/a£=Ky,  P^/a?  and  KTTy  as  functions  of  r are  shown  in  Figure  4. 
It  is  seen  chjjt  the  normalized  timing  error  variance  just  after  processing 
a data  point  P i/af  is  always  less  than  1,  while  the  normalized  timing 
error  variance  just  prior  to  processing  a data  point  p-|/a?  can  be  larger 
or  smaller  than  1.  The  cross  over  point  occurs  at  r % 16.  Therefore  if 
the  observation  interval  is  chosen  to  be 
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(74) 


Tf  " • 


then  the  timing  error  can  be  kept  below  the  inherent  measurement  error 
at  all  times. 

A simplified  version  of  the  SDDLL  including  a discrete  Kalman  filter 
is  shown  in  Figure  5.  Referring  to  Figures  3 and  5,  k]  is  the  pro- 


e 


Figure  5.  Model  of  the  SDDLL  with  a discrete  Kalman  filter. 


portionality  constant  of  the  sampled  error  voltage  versus  timing  error 
when  | £ | < a/2,  i .e. , 

Es  = k-|£  . (75) 

The  time  delay  e"^Ds  accounts  for  the  finite  response  time  of  the  band- 
pass filter  in  the  SDDLL.  The  integrator  k£/s  is  the  digital  voltage 
controlled  clock.  The  sampler  Tf/N  is  provided  because  time  corrections 
can  only  be  made  discretely.  The  value  of  N should  be  large  enough 
that  timing  error  due  to  the  range  rate  will  be  a small  fraction  of 
the  chip  length  and  the  clock  will  remain  locked. 

The  samplers  labeled  Tf  in  Figure  5 close  momentarily  once  per  sample 
period  at  times  t = nTf+Tg  where  nTf  indicates  the  arrival  times  of  the 
samples  T,  while  the  sampler  labeled  Tf/N  closes  N times  per  sample 
period  (evenly  spaced)  with  the  initial  closure  occurring  at  t = nTf  + Tq 
+ Tf/2n. 

The  filter  output,  Tc  in  Figure  5,  provides  an  optimum  timing 
estimate  (within  limits  established  by  the  discrete  corrections 
Jf/N)  at  all  times.  This  estimate  varies  from  the  predicted  estimate 
Tn  (with  variance  Pn  - the  worst  case)  just  prior  to  processing  the 
nth  data  point  to  the  filtered  estimate  Tn  (with  variance  Pn  - the 
best  case)  just  after  the  nth  data  point  has  been  processed. 
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If  the  bandwidth  of  the  band-pass  filter  is  optimum  [1],  i.e.. 


Btf  = 1.37/Ma 
^ 'opt 

where  Ma  is  the  length  of  the  PN  code,  then 

|l/2 


. 1 

1“  " ? 


1.215  . 1.476 


where 


)2  for  square-law  detection 
|l  for  linear  detection 


E = P^Ma  = input  pulse  energy  (joules) 


Pr  = input  power  (watts) 


and 


(76) 


(77) 


(78) 


Nq  = one-sided  input  noise  power  spectral  density  (watts  per  Hertz). 

To  obtain  the  standard  deviation  of  the  maneuver  noise  oa,  we 
assume  that  the  probability  density  of  T is  as  shown  in  Figure  6 [10]. 


O 

CL 

p 

P 

i m 

i m 

► 

-at  o at 

Figure  6.  Probability  density  of  T. 

The  quantity  Ay  in  this  figure  might  be  considered  as  the  maximum 
"path  delay  acceleration"  given  by 
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(79) 


At  = A/c  (sec/sec; 


where 


p 

A = terminal  acceleration  (m/sec  ) relative  to  the  satellite,  and 
c = velocity  of  light  (m/sec). 

Therefore,  since  the  mean  value  E(T)  is  zero  by  symmetry,  the  variance 
is  given  by 

°°  A2 

°a  = E(Y2)  = J iZ  p(Y)  O + 4pm-P0)  (sec"2)  (80) 

_ oc 

where  p(T)  is  the  probability  density  function  of  Figure  6.  The  value  of 
oa  thus  obtained  is  taken  to  be  a known  constant  in  the  modeling  of  the 
Kalman  filter. 

If  or  or  aa  vary  widely  with  time  in  a given  application,  separate 
measurements  (or  estimates)  of  o^/oa  could  be  made  periodically  and  used 
in  place  of  the  calculated  constants  to  obtain  improved  filter  performance. 

C.  Augmented  Discrete  Kalman  Filter 

In  the  previous  section  the  maneuver  noise  T is  considered  to 
be  white.  However,  in  a TDMA  modem  a typical  sampling  period  Tf 
is  only  about  ,01  sec  which  is  much  shorter  than  the  time  required 
for  executing  a manuever,  hence  the  maneuver  at  one  sampling 
period  will  be  correlated  with  the  maneuver  at  the  previous  and  at  the 
next  sampling  period.  In  fact  it  can  be  correlated  with  maneuvers 
separated  by  many  sampling  periods.  If  such  is  the  case,  the  simple 
discrete  Kalman  filter  discussed  in  the  previous  section  may  not  yield 
satisfactory  performance.  An  augmented  Kalman  filter  is  therefore 
needed. 

First,  Equations  (12)-(14)  must  be  modified.  A satisfactory  cor- 
relation function  in  the  sense  of  providing  realistic  approximations 
to  the  statistics  of  the  maneuvers  as  well  as  yielding  mathematically 
tractable  solutions  is  [10] 

E[T(n)  T(k)]  = a2p'n'k'  (81) 

where  at  is  obtained  as  before  (Figure  6 and  Equation  (36))  and  the 
correlation  coefficient  p is  modeled  by 
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1 - AT 


if 


Tf  ix 


if  W 


(82) 


The  quantity  A is  essentially  the  inverse  of  the  average  manuever  duration. 
With  the  correlation  defined  as  above,  the  dynamics  of  the  terminal  are 
then  described  by 


T(n+1 ) = T(n)  + f(n)  Tf  + J-  T(n)  T2f 
t(n+l)  = f (n)  + T(n)  Tf 
T(n+1 ) = p T (n)  + u(n) 

2 2 

where  u(n)  is  a white  noise  of  zero-mean  and  variance  aa(l-p  ).  In 
matrix  notation,  ( 83 ) - ( 85 ) becomes 


(83) 

(84) 

(85) 


Z(n+1 ) = «Z(n)  + ru(n) 


where 


Z(n)  = 


4 = 


and 


r = 


T(n ) 
T(n) 
T(n) 


1 Tf  T-/2 

0 1 T 

0 0 p 


(86) 


(87) 


(88) 


The  measurement  is  given  by 
X (n)  = H Z(n)  + c(n) 

wi  th 

H = [1  0 0] 

and  the  statistics  of  s(n)  as  given  previously  in  Equations  (7)-(9). 


(89) 

(90) 

(91) 
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Mathematical  development  of  the  augmented  Kalman  filter  fnllnwc 
essentially  the  same  path  as  in  the  previous  section  a f°"°WS 
changes  in  the  matrix  contents  (see  Equations  (28)  through  (36)), 


-fO-P2) 


Manipulations  for  solving  this  system  of  equations  are  as  follows: 


(99)  - | x (98): 

D = | B2/(l+A-p)  . 


Substituting  (109)  into  (101): 

F = . 4p2B4 

r2( 1+A-p )2( 1+A) ( 1-p2) 

Substituting  (110)  into  (98): 

1 B2  2p2B4 

2 2(1+A)  r2( 1+A-p )2(1+A)( 1-p2) 

Substituting  (109),  (110)  and  (111)  into  (100): 


B 


2 


2pB 

r(l+7r-p)T1-p}_ 


2 


(1+A) 


or 


1 + 


2PB 

r(l +A-p )( 1-p  ) 


(96)  + i x (97)  + x (98): 

r 


A - i B = (A  + i B - B2)/(l+A) 

or 

B = J (/T+A  - 1 )2  . 


(109) 


(no) 


OH) 


(112) 


(113) 


(114) 


(115) 


28 


Now  let 


and 


A = (x+ir-1 
r = r (^)1/2 

p 1-p' 


(116) 

(117) 


Then  substituting  (116),  (116)  and  (117)  into  (113)  gives 


x4+2(l-p )(+l-  —■) x3+( 1 -p  ) 

p 


(1-p)*  y~ 
pj 


x2_  2( 1-p )(3-P ) x+  2(1 -p)d  = 

r r 

p p 

(118) 


It  can  be  shown  (since  A = ^/ar  must  be  real  and  non-negative)  that  the 
lower  sign  options  of  Equation  f 1 1 8 ) must  be  used,  and  in  addition  that 
this  equation  yields  only  one  positive  real  root,  which  is  the  correct 
root. 


The  4th  order  equation  (118)  may  be  solved  by  the  standard  (closed 
form)  method  as  given  in  many  mathematical  handbooks  [11-12]  or  by  computer 
iteration  techniques.  Both  methods  have  been  used  to  generate  results 
presented  in  this  report.  (The  iteration  techniques  are  easier  to 
program  while  the  closed  form  solutions  are  more  efficient.) 


o 

If  r^(l-p)  >>1,  then  an  approximate  solution  for  x is 


(119) 


From  (34)  and  (38)  we  obtain  the  steady-state  gain  matrix. 


kt 

2 

P1 

2 

Kf 

= 

p2 

- . 2 

_kt_ 

where 


(120) 
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Kt  - P,/oJ  - 


Kf  = P2/oJ 


-v,  ? 

1 + P,/o‘ 


pp 


V..  _ n / 2 _ 4p 

Kt  ■ V° 5 ■ ~3 

rTf 


fy0e0aTf 
~ — 2 — 
1+V°5 

Vp<Va 

^VT 


(121) 


(122) 


(123) 


Then  we  obtain  from  Equations  (96)  to  (123) 

-v, 


A = (x+1 )2-l  = 

(124) 

(125) 

? a f 

2 p 

n - 2B  - 4 

r(l+A-p ) poo 

t,  a 

(126) 

K ^ - A 

T " T “ 1+A 

(127) 

X.T  - 4 8 

KTTf  r T+ft 

(128) 

u..t2  _ 4p  D 

Vf ' ~ we  • 

(129) 

Graphs  of  frj/o2,  P-j/o2  = Ky,  K|Tf, 

2 

and  twTr  as  functions  of  r 

T f p 

for  different  values  of  p (obtained  from  the  exact  solution.  Equation 
(118))  are  shown  in  Figures  7 through  10.  Note  that  for  p=0  the  augmented 
Kalman  filter  reduces  to  the  two-dimensional  filter  analyzed  earlier,  and 
the  corresponding  curves  are  in  agreement  wi th  those  given  previously  in 
Figure  4. 
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A model  of  the  resulting  augmented  Kalman  filter  incorporated 
into  the  receive  clock  sampled  data  delay  lock  loop  is  given  in 
Figure  11.  Sampling  switch  timing  for  this  model  and  comments  con 


° — ►/eY-J  k.e~Tos 


I kt /k. k. 


kt/  k|  k2 


Ky/k.kj 


Figure  11.  Model  of  the  SDDLL  with  an  augmented  Kalman  filter. 

cerning  the  time  delay  and  the  output  Tc  are  the  same  as  that  previously 
given  for  the  tv/o -dimensional  model  of  Figure  5. 

D.  Transient  Response  of  the  Discrete  Kalman  Filter 


The  settling  time  is  an  important  factor  to  consider  when  choosing 
a suitable  filter  for  implementation  in  the  SDDLL.  A Kalman  filter, 
though  optimum  in  a sense,  requires  more  computer  storage  and  computer 
time  than,  for  example,  the  (steady  state)  Wiener  filter.  If  the  settling 
time  is  sufficiently  short,  a corresponding  Wiener  filter  can  be  used  in 
place  of  the  more  sophisticated  Kalman  filter.  In  this  section,  the 
transient  responses  of  the  discrete  Kalman  filters  discussed  in  previous 
sections  will  be  studied  and  curves  of  the  transient  response  and 
approximate  expressions  for  the  settling  time  will  be  obtained. 


The  transient  response  of  the  simple  discrete  Kalman  filter  is 
studied  first.  Combining  Equations  (28)  and  (29), 

• _1[?((n+l  )-rqrT]*'1T  = fan)  - ^'(n)HT[HP(n)HT+R]"1H  ^(n). 


(130) 


Equating  corresponding  elements  on  both  sides  of  (130),  we  obtain 


P-,  (n+1 ) -2^(n+l)Tf+^3(n+l)T^-0^/4  = ^ (n)o^/  (n)+a*)  (131) 

£2(n+l)  - P3(n+l)Tf  + a^/2  = ^(n)o^/ (^  (n )+o^)  (132) 


£3(n+l)  - a2J2f  = P3(n)  - ^(n)/(^(n)+o^)  (133) 

'V  X % % 

where  Pq , ?2,  P3  are  the  elements  of  P as  defined  previously  in  Equation 
(40).  In  the  steady  state,  Equations  ( 1 31 ) - ( 1 33 ) reduce  to  Equations 
(48)-(50).  The  above  difference  equations  are  nonlinear  and  can  not  be 
solved  analytically,  hence-  numerical  methods  must  be  employed.  To  start 
the  iterations,  we  obtain  the  initial  value  of  P(n)  from  (28)  and  (36) 


$(3)  = 


, ? ? 4 

5(a^+a^/10) 


<°KT4f/  4) 


3 / 2 2,4..,  2 , 2.,  2T4 /Q, 

T-  (^+aJf/4)  ^ (a5+5oaTf/8) 

It  is  convenient  to  define  the  following  variables; 

r=V°aTf 

A(n)=?.j  (n)/o^ 

B(n)=P2(n)/a^oaTf 


C(n)=?3(n)/a^. 


Then  (131)  through  (134)  become 

A(n+1)  = B ( n+1 ) - C ( n+1 ) + + 

B(n+1)  = i C(n+1 ) + | 


(134) 


035) 

(136) 

(137) 

(138) 

(139) 

(140) 
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* c<">  - -®)  * 1 


(141) 


C • h)< 

r(!  * ?)°<Vf 


r 1 + 


1 21,  . 10' 
ff  r ( 1 1 * 1 

K r . 


> 2T2 
j °aTf 


(142) 


The  filtered  covariance  elements  are  related  to  A(n),  B(n)  and  C(n) 
by 

P-j  (n)/o^  = A ( n ) / ( 1 +A ( n ) ) 

P2(n)/a5oaTf  = B(n)/(1+A(n) ) 

^3(n)/aaTf  = c(")  ' B2(n)/(1+A(n) ) 
and  the  gain  factors  are 


(143) 


(144) 


(145) 


KT(n)  = P^nJ/oj 


(146) 


Kt(n)  = rT^  ' ~P2(n)/VaTf* 


(147) 


* 2 

Graphs  of  P-|(n)/a^,  K^.(n)  and  K|(n)T^  for  different  values  of  r are 

given  by  the  e=0  curves  in  Figures  12,  13  and  14.  From  these  figures, 
we  obtain  a fairly  good  approximation  for  the  settling  time, 


At  % 2.5  r,Z  T 


(148) 
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Figure  13.  Normalized  filtered  estimate  error  variance  and  Kalman  gain 
transient  behavior  for  the  Kalman  filter. 


For  example,  using  a typical  value  of  Tf  % .01  sec 


At  K, 


.08  sec  , 
/.  8 sec  , 
8 sec  , 


r=10 

r=103 * 

r=105 


0 49) 


From  the  definition  of  r (Equation  (135)),  Equation  (148)  indicates 
that  the  settling  time  depends  only  on  the  ratio  or/oa,  irrespective  of 
the  sampling  period  Tf,  However,  we  have  neglected  correlation  in  the 
maneuver  noise  when  deriving  (148).  As  the  sampling  period  becomes  shorter, 
the  effect  of  correlation  in  the  maneuver  noise  can  not  be  overlooked.  We 
expect  that  the  settling  time  will  also  depend  on  the  sampling  period 
when  correlation  is  included. 


Now  consider  the  transient  response  of  the  augmented  discrete 
Kalman  filter.  Following  the  same  reasoning  as  above,  we  obtain  the 
recursive  covariance  equations  as  follows: 


A(n+1 ) = B(n+1 ) - 


16 

T 

r 


C(n+1 ) - £ D(n+1 ) 


16 

T 


E( n+1 ) - ^ F(n+1) 


+ 


A(n)  4 
T+ATriy  ■ -7 


(150) 


B(n+1 ) 


7 C(n+1 ) + D( n+1 ) 


- “ E( n+1 ) + - F ( n+1 ) 


B(n 
1 +A ( n ) 


+ 


2 

r 


(151) 


2 

C (n+1 ) = 2E(n+l)-F(n+l)+C(n)  - - 1 (152) 

D(n+1)  = i E(n+1)  - £ F(n+1)  + \ (153) 

E(n+1 ) = F ( n+1 ) + PE(n)  - p^|°jn)  + 1 (154) 

2 2 

F( n+1 ) = P2F(n)  - (155) 
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where  r,  A(n),  B(n)  and  C(n)  are  defined  as  before  and 


D(n)  = P^(n)/po^aa 


(156) 


E(n)  = P5(n)/po^Tf 


F(n)  = (^6(n)-Og)/p2o^ 

In  the  steady  state  ( 1 50 ) - ( 1 5 5 ) reduce  to  (96)-(101). 
variance  matrix  is  obtained  from  (94)  and  (28) 


(157) 


(158) 


The  initial  co- 


P(3)  = 


5 „ + MJM)  2 

5r  / 5 


3 + 4(1^P  An  a T 2(1 +p  ) -, 

V r2  / C aTf  r~^  C < 


3 r + iil±dcaaTf  1 r2  (l  + cjft  (l+  U pa2Tf 


2 ( 1 +p  ) 

r^Va 


(1+  ^ P)  POgTf 


(159) 


The  filtered  covariance  matrix  elements  are  then  given  by 


P-,  (n)/a2  = A ( n ) / ( 1 +A ( n ) ) 

(160) 

P2(n)/a?aaTf  = B(n)/(1+A(n  ) ) 

(161) 

P3(n)/o2T2  = C(n)-B^n)/(1+A(n)) 

(162) 

P4(n)/po_aa  = D(n)/(1+A(n) ) 

(163) 

P(:(n)/pa?Tf  = E(n)  - B(n)D(n)/(l+A(n)) 

0 a T 

(164) 

P6(n)/p2aa  = F(n)  - D2(n)/(1+A(n))  + 1/p2 

(165) 
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and  the  gain  factors  are 


K-j-(n)  = P1  (n  )/a2  (166) 

Kf  (n)  = P2(n)/a^oaTf  O67) 

Kj ( n ) = — ^ P4(n)/po  o . (168) 

rTf 

x 2 2 

Graphs  of  Pj (n)/af  , K-j-(n),  K|(n)T^  and  Ky(n)T^  for  different  values 

of  r and  p are  shown  in  Figures  12  through  15.  From  these  figures  we 
obtain  an  approximation  for  the  settling  time 


At  * 2.5r1/2  Tf  % 2.5  r.1/2Tf,  if  r1/2(l-p)>l.  (169) 

The  condition  r^ 72( 1 -p  )> 1 is  essentially  the  same  one  under  which  (119) 
is  obtained.  For  fixed  r,  the  settling  time  decreases  as  p increases 
until  the  condition  rl/2(l-p  )--l  is  no  longer  satisfied.  Then  the 
settling  time  increases  as  p increases  further.  The  dependence  of  At  on 
the  sampling  period  is  implicitly  contained  in  p (p=l-ATf).  For  r<10^, 
p<.999  and  Tf = . 01  sec,  the  settling  time  falls  approximately  in  the  range 
of  one-tenth  to  ten  seconds. 


E.  Monte  Carlo  Simulation 


Performance  of  the  discrete  Kalman  filter  described  in  the  previous 
sections  can  be  investigated  using  a Monte  Carlo  simulation  method.  In 
this  method,  the  receiver  noise  and  the  maneuver  noise  are  generated  by  a 
computer  using  arithmetical  random  number  formulas.  The  results  of 
several  simulation  studies  are  summarized  here. 

In  the  first  case,  the  simulation  method  was  used  to  obtain  a 
comparison  between  the  performances  of  the  discrete  Kalman  filter  and 
its  corresponding  Wiener  filter.  Correct  statistics  of  the  receiver 
noise  and  the  maneuver  noise  were  assumed.  Values  of  p,  oa  and  o£  were 
chosen  as 
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r the  Kalman  filter 


p = 0,  .5,  .9,  .99,  .999 

O70) 

o = .5A/sec^,  15A/sec^ 

Q 

(171) 

and 

0 = .01  A,  .1  A. 

(172) 

2 

For  a code  chip  rate  of  40  Mbit/s,  a = 25  nsec.  Then  aa  = 15A/sec 
represents  an  acceleration  of  llg  which  is  probably  the  highest  for 
aircraft.  For  oa  equal  to  15A/sec^,  the  sampling  period  Tf  should  not 
be  greater  than  .02  sec.  in  order  to  keep  the  transient  timing  error 
below  0.5a  and  thus  avoid  loss  of  lock.  In  this  simulation  study  Tf 
was  chosen  to  be  .01  sec. 

In  the  steady  state,  the  standard  deviation  of  the  one-step- 

prediction  timing  error  was  found  to  agree  in  general  to  within  10% 
of  the  analytical  results  of  Section  C,  except  for  aa  = .5A/sec2, 
o ^ = .1a,  in  which  case  the  simulation  results  were  30%  lower  than  the 
analytical  results. 

The  settling  time  of  the  discrete  Kalman  filter  agrees  with 
Equation  (169).  For  the  corresponding  Wiener  filter,  the  settling  time 
is  about- twice  as  long.  For  =.1a  the  timing  error  of  the  Wiener 
filter  in  the  transient  state  can  become  greater  than  .5a  which  could 
unlock  the  local  clock  of  the  SDDLL.  This  is  also  true  for  the  Kalman 
filter  if  > .2a.  It  appears  that  a better  initial  estimate  than  that 
given  by  Equations  (93)  and  (94)  is  needed; if  a^,>  ,2a.  The  crucial 
points  are  the  initial  estimates  of  T(n)  and  of  T(n).  Error  in  the 
initial  estimate  of  T(n)  is  rather  unimportant  because  of  the  short 
sampling  period. 

In  the  second  simulation  study,  the  performance  of  the  simple 
discrete  Kalman  filter  was  compared  to  that  of  the  augmented  Kalman 
filter  of  the  previous  study.  If  a simple  discrete  Kalman  filter  is 
implemented  irrespective  of  the  correlation  coefficient,  the  per- 
formance will  not  be  optimum.  The  higher  the  value  of  p,  the  poorer 
the  performance  will  be. 

Assuming  correct  values  of  cr,  and  Jr  in  implementing  the  simple 
Kalman  filter,  the  performance  was  degraded  by  3%  for  p = .5,  40%  for 
p = .9  and  60%  for  p = .99.  However,  the  analysis  of  previous  sections 
indicates  tu’t  the  maneuver  noise  standard  deviation  oa  should  be  re- 
placed by  aaL(l+p)/(l-p)]V2  if  a simple  discrete  Kalman  filter  is  to  be 
implemented  for  p/0  (i.e.,  this  is  equivalent  to  replacing  the  parameter 
r by  rr  - see  Equations  (102)  and  (117).  Improvements  in  performance 
can  then  be  obtained.  This  was  verified  in  the  simulation  study  except 
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for  the  case  p = .99,  oa  = 15A/sec^,  = .01A.  For  other  parametric 

values  the  degradation  in  the  performance  was  only  1%  for  p = .5,  10% 
for  p = .9  and  30%  for  p = .99.  Thus  it  appears  that  with  p less  than 
.99  and  with  proper  substitution  of  the  maneuver  noise  variance,  the 
simple  discrete  Kalman  filter  can  be  implemented  to  save  computation 
time  and  storage  cost  while  still  providing  satisfactory  performance. 

For  p equal  to  .99  and  higher,  the  augmented  discrete  Kalman  filter  should 
be  used. 

One  potential  major  disadvantage  of  implementing  the  discrete  Kalman 
filter  is  that  it  may  require  accurate  estimation  of  the  value 
r = 4o^/oaT^.  Therefore  in  the  third  simulation  study  the  effect  of  the 
accuracy  of  r on  the  performance  of  the  discrete  Kalman  filter  was 
studied.  First,  for  p = .99  and  .999,  it  was  found  that  the  performance 
of  the  augmented  filter  is  degraded  by  only  10%  if  the  estimation  of  r 
is  off  by  a factor  of  2,  and  by  80%  if  the  estimation'  of  r is  off  by  a 
factor  of  5.  Testing  values  of  oa  and  are  the  same  as  in  the  first 
two  studies.  The  small  percentage  of  deg'radation  can  be  attributed  to 
the  small  slopes  of  the  curves  shown  in  Figure  7.  The  same  study  was  then 
repeated  for  p = 0,  .5  and  .9  using  the  simple  discrete  Kalman  filter  with 
proper  substitution  of  the  manuever  noise.  Results  were  essentially  the 
same  as  for  the  augmented  filter  with  only  a slight  increase  in  the 
percentage  degradation  of  performance. 

The  tracking  capability  of  the  discrete  Kalman  filter  (worst  case) 
was  then  considered.  The  terminal  was  initially  stationary  and  the 
discrete  Kalman  filter  with  assumed  values  of  p,  og  and  a g was  in  the 
steady  state.  Then  an  instantaneous  velocity  was  suddenly  imparted  to  the 
terminal  as  when  a frequency  offset  occurs.  Typical  response  of  the 
discrete  Kalman  filter  is  shown  in  Figure  16.  Following  an  initial 
over-shoot  the  filter  settles  down  quickly.  It  was  found  that  the 
maximum  instantaneous  velocity  which  can  be  imparted  to  the  terminal 
without  risk  of  unlocking  the' local  clock  of  the  SDDLL  is  approximately 
given  by 

’-80  r-1/2  VTf  . (173) 

Results  of  additional  simulation  studies  of  the  discrete  Kalman 
filter,  including  transient  response,  noise  performance  for  stationary 
terminals,  and  comparisons  of  analytical  values  of  error  variance  with 
simulation  results  are  given  in  Section  VII. 

F.  Up-Link  Synchronization 

Closed-loop  synchronization  of  the  TDMA  system  is  achieved  by 
adjusting  the  timing  of  up-link  transmissions  to  maintain  a prescribed 
timing  relationship  between  the  received  down  link  signals  and  a 
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network  clock  signal.  The  synchronization  of  the  single  one  way  link 
illustrated  in  Figure  17  [1]  will  be  considered  in  this  section. 


O SATELLITE 


cp 

CO 


0 - DATA  TRANSMITTING  (ORIGINATING)  TERMINAL 
D - DATA  RECEIVING  (DESTINATION)  TERMINAL 


Figure  17.  A single  (one-way)  TDMA  link. 

As  previously  stated,  the  terminals  need  not  be  ground-based.  Fast 
moving  terminals  can  be  included.  Following  [1],  it  will  be  assumed 
that  a portion  of  the  time  allocated  to  each  terminal  is  devoted 
exclusively  to  transmission  of  a pulsed  envelope  ranging  signal  from 
which  timing  error  information  is  extracted.  Strict  preservation  of 
the  TDMA  concept  requires  that  the  network  clock  signal  also  have  a 
pulsed  envelope  waveform. 
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The  TDMA  synchronizing  receiver  employs  two  sampled-data  delay- 
lock  loops  (SDDLL's)  to  maintain  proper  signal  timing.  In  one  SDDLL, 
the  time  base  of  a locally  generated  clock  signal  is  aligned  with  the 
time  base  of  the  received  network-clock  signal.  This  SDDLL  with  a 
discrete  Kalman  filter  has  been  analyzed  in  previous  sections.  The 
second  SDDLL  estimates  the  timing  error  of  the  received  ranging  signal 
relative  to  the  locked  local  clock  and  appropriately  corrects  the 
timing  of  the  up-link  transmissions.  The  performance  of  the  second 
SDDLL  will  be  discussed  in  this  section.  Both  the  network  clock  and 
ranging  signal  waveforms  are  bi-phase  modulated  by  suitable  PN  codes. 

An  abbreviated  block  diagram  of  the  synchronizing  receiver  is  shown  in 
Figure  18  [1],  (Identical  codes  for  the  clock  and  ranging  loops  are 
assumed  here  for  simplicity,  although  different  codes  are  actually  used 
for  signal  identification.)  Note  that  the  averaging  loop  filter  in  the 
clock  loop  of  the  original  TDMA  modem  has  been  replaced  by  a discrete 
Kalman  filter.  However,  as  stated  previously  in  Section  II,  within  the 
spectral  range  of  the  timing  signal  and  the  dynamic  range  of  the  terminal 
motion,  the  discrete  Kalman  filter  in  the  clock  loop  can  provide  the 
tracking  function  for  both  loops.  Therefore,  in  the  ranging  loop  the 
simple  averaging  loop  filter  will  be  retained. 

Closed-loop  timing  corrections  in  the  ranging  loop  must  be  separated 
by  an  interval  larger  than  the  propagation  delay  of  the  ranging  signal 
(about  250  msec  when  a synchronous  satellite  relay  is  used)  in  order  to 
maintain  ranging  loop  stability.  However,  if  fast  moving  terminals  and 
a high  code  chip  rate  are  employed,  timing  error  in  the  ranging  loop 
accumulated  in  such  a time  interval  can  be  larger  than  A/2  if  only 
closed-loop  corrections  are  made.  Therefore,  estimates  of  the  down-link 
delay  rate  and  acceleration  derived  in  the  clock  loop  are  employed  as 
open-loop  corrections  in  the  ranging  loop.  The  modified  two-loop  TDMA 
synchronizer  used  previously  [1],  with  certain  changes  as  shown  in 
Figure  19,  will  serve  this  purpose.  In  this  figure,  ATcf  and  ATrf  are 
the  sampling  period  offsets  (errors)  in  the  clock  loop  and  the  ranging 
loop,  respectively  (i.e.,  ATc>rf  = Tc>rf-Tf  where  Tc  rf  are  the  nominal 
subframe  lengths  and  Tf  is  the  actual ’subframe  length  at  the  satellite). 
The  formula  for  determining  the  quantity  A0  will  be  given  later.  Details 
of  the  clock  loop  Kalman  filter  were  given  previously  in  Figure  11. 

Relationships  between  the  sampling  instants  of  the  samplers  are 
shown  in  Figure  20,  where  Ta  is  the  PN  code  period.  It  will  be  assumed 
that  the  Nr  ranging  signal  pulses  are  sampled  every  Tf  seconds  where 
Tf  = KTa.  Neglecting  circuit  delays,  the  sampling  instants  (A)  represent 
the  times  when  system  clock  samples  are  received  and  corrections  in  local 
clock  timing  are  made  as  a result  of  measured  timing  errors.  At  times 
(B)  additional  corrections  are  made,  as  predicted  by  the  Kalman  filter, 
to  compensate  for  the  effects  of  terminal  or  satellite  maneuvers 
(velocity  and  acceleration).  Times  (C)  indicate  when  ranging  pulses  are 
received  on  the  down-link.  Note  that  the  Nr  ranging  pulses  of  each  burst 
are  spaced  by  the  subframe  time,  Tf,  as  are  the  system  clock  pulses  but 


51 


SATELLITE 


CLOCK  LOOP  COOE  GENERATOR 


Figure  20.  Sampling  instants  for  the  case  Nr=2,  and  Nc=2. 

are  offset  from  these  by  an  amount  nr  Ta  (when  properly  timed)  so  that 
interference  between  system  clock  and  ranging  signals  does  not  occur. 

Note  also  that  to  prevent  instability  in  the  ranging  loop,  Tf  and  Ir 
(the  number  of  subframes  per  frame)  are  chosen  so  that  the  separation 
between  the  first  ranging  pulse  of  each  burst  and  the  last  pulse  of  the 
previous  burst  exceeds  the  maximum  round-trip  propagation  delay,  Tmax. 
Sampling  instants  (D)  represent  the  times  when  the  Nr  ranging  pulses  of 
each  burst  have  been  averaged  and  corrections  in  transmit  timing  to 
compensate  for  measured  timing  errors  are  made.  Additional  (open-loop) 
corrections  to  transmit  timing  are  also  made  at  times  (A)  and  (B) 
above  (see  Figure  19). 

We  now  proceed  to  analyze  the  responses  of  the  modified  synchronizer 
as  follows.  Within  a sampling  interval  Tf,  the  network  clock  can  be 
described  by 


To(nTf+6)  = T0(nTf)  + fQ(nTf)6  + T0(nTf)«V2,  0<S<Jf. 


(174) 


In  the  steady-state,  the  discrete  Kalman  filter  can  track  non-random 
delay  accurately.  Therefore,  in  the  steady-state 
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(175) 


E[Tc(nTf+5)]  = E[To(nTf)]  + E[tQ(nTf)]6  + E[TQ(nTf ) ]«2/2. 

where  E[  ] indicates  the  expected  value.  However,  this  equation  is 
true  only  when  the  corrections  from  samplers  (B)  of  Figure  19  can  be 
made  continuously.  To  account  for  the  discrete  nature  of  these 
corrections,  (175)  is  modified: 

Tf  . Tf 

E[Tc(nTf+m  p-  +e)]  = E[T0(nTf)]  + E[t  (nTf)]  m p- 
c c 

? Tf  ATrf 

+ E[T  (nTf )]ni  ^ ^ 6. 

° t 2rr  f 

C 

where 


(176) 


Tf  Tf  Tf 

- ^j—  < e < jip  > 0 1 m < Nc  and  *-)  <_  m p-  + 6 T ^ . (177) 

c c c 

Thus  the  expected  value  of  the  clock-loop  timing  error  is 

Tf  Tf  Tf 

ecs^nTf+m  n“  + = ECVnV" } r +e^  ■ Tc(nTf+m  1 — + 

c c c . 

-|E[t0(nIf)]  - ^ * E[T0(nTf)]  (\  m ^>e  . 

J (178) 

The  subscript  s denotes  the  expected  value.  From  the  previous  sections 
we  obtain  the  variance  of  the  clock-loop  timing  error, 

a\  (nTf+6 ) = + 2P2«  + P3 62  + P^2  + P^3  + F>664/4,  6 = m p- 


and 


\ (nTf+Tf-0)=?f1 


(179) 


(180) 
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where  the  P^'s  and  l^'s  are  elements  of  the  covariance  matrices  P 
and  P,  respectively,  as  defined  in  (95).  The  variance  c>2  (nTf+6)  is 
a monotonical ly  increasing  function  of  6 . Ec 

Next,  consider  the  ranging-loop  response.  We  shall  first  derive 
some  additional  formulas  related  to  the  steady-state  Kalman  filter 
which  are  essential  to  the  calculation  of  the  ranging-loop  timing 
error.  From  (86), 

SL 

Z(n+£)  = ? Z(n)  + £ «*"k  ru(n+-l)  (181) 

k=l 


where 


1 T, 


* = 0 1 


0 0 


(182) 


i n 

1 iTf  T [Pk_1(  2*-2k+l)  1 Uf  /[2Ml-o)-(l+p)(l-pi)/(l-p)2] 


• = 0 1 


0 0 


T V _ k— 1 

f khp 


■ 0 1 Tf(l-e  )/(1-p) 


0 0 


(183) 


TJn+i)  = pl  T (n)  + [ p*'k  u(n+k-l) 
0 0 k=l 


(184) 


Vn+‘)  = Vn)  + Tf  [j^T  V")  + j,  T^p~~  u(n+k-]) 


(185) 


"" 


T0(n+i) 


T0wmf  t0(n)  * ^M’-al^£luV 

U_P  ) 


V"> 


L 

k=l 


2(^-k)(1-p)-(Hp)d-P^k) 


(1-P)‘ 


(186) 


• u(n+k-l)| 

Also  from  (19),  (21),  (22)  and  (86), 


Z (n+1)  = Z(n+1  )-Z(n+l  ) = ( I-KJ-I)  [i=Z  (n)+l'u(n)  ]-K€  (n+1 ) 


(187) 


and 


Z (n+*)  = [ (I -KH)®]XZ  (n)  + \ [ ( I - KH ) * [ ( I - KH ) ru ( n+k- 1 )-K?  (n+k)] 


e'  ' k=l"  ' c 


(for  l > 0).  (188) 


Thus 


P(*)  2 E[Z£(n+Jl)Z^(n)]  = [(I-KH)*]AP 

oe(A)  = E[Z£(n+£)5c(n)]  = -[( I-KH)  *]V2 

oU)  = E[Z  (n«)u(n)]  = [( I-KH)*]*"1  ( I-KH)  Ta^(  l-P2) , 

a £ ci 


(189) 

(190) 

(191) 


Note  that  a r{SL)  and  o («■)  are  vector  quantities,  i.e., 
5 a 


o U)  = 

a 


0aT(i) 

°af(i) 

aaT(t) 


(192) 


and  similarly,  for  a^U). 
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Now  refer  to  Figure  19.  The  ranging-loop  clock  at  IrTf+£Tf+nrTa 
(i.e.,  the  time  at  which  a ranging  pulse  is  received)  is  related  to  that 
same  clock  at  £Tf+nrTa  (one  frame  earlier)  by  the  expression 


AT 


T (I  T,+«.T,+n  T ) = T (*T,+n  T ) 

K'  rT  7"  Pel  p I r a' 


rf 

t: 


I T, 
r f 


+ C U,I  +*)  + C (a, I +1)  + C (193) 

t r o r c 

where  CT,  CQ  and  C are  due  to  the  propagation  delay  change,  the  open- 
loop  correction  ana  the  closed-loop  correction,  respecti vely. 

Because  of  the  propagation  delay,  any  change  in  satellite  motion  is 
not  known  at  the  terminal  until  after  the  (one-way)  propagation  delay 
has  elapsed,  while  a change  in  terminal  motion  is  known  immediately. 

Hence  the  effects  of  the  propagation  delay  change  CT  differ  depending 
upon  whether  the  satellite  or  the  terminal  is  maneuvering.  The  case 
of  a stationary  terminal  and  maneuvering  satellite  (worst  case)  will  be 
considered  first.  In  this  case  CT  is  given  by 


CT(Mr+A)  = 2[T0(IrTf+*Tf+nrTa) 


VWa^ 


(194) 


Expressions  for  the  open-loop  and  the  closed-loop  corrections  involve 
the  propagation  delay  t and  the  ranging-loop  sampling  instant  nrTa. 


Let 


T 


n T 
r a 


N T 

T f 


e 

T 


(195) 


i.e.,  the  time  period  from  transmission  of  the  ranging  pulse  to  the 
nearest  subframe  boundary  (and  received  system  clock  pulse)  preceeding 

reception  of  this  pulse  is  expressed  as  an  integer  number  (NT)  of  sub- 
frames  (Tf)  minus  an  integer  number  (mT)  of  clock  loop  correction  inter- 
vals (Tf/Nc)  minus  the  remaining  offset  (eT), 


where 


Tf  Tf  Tf 

" 21T  < 6t  < ’ 0 1 mT  - Nc  and  0 - mT  W + \ - Tf*  (196) 

c c c 
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C0(Mr+*>  - - [Tc(IrTfmf-Vf)-Tc(*Tf-NJf)  - IrTfl 


/III  I r 

■ttc(IrTf+tV",Tf>-tc(<TrN,Tf)H-frI  + flo 


fm  Tf  "\  m Tf 

- [Tc(,r1ft‘Tf-'l,'If>-Tc(1TrN,Tf)^5fi7  + V ^7 


(197) 


A.  V 


Cc  ‘ if  L <Tc(nTftnrTa>-«r<nWa)-Tr<nTftnrTa) 
r n=0 


(198) 


x+iNr-l  )Tf-IrTf  <*  Tf  < t + (Nr-DTf 


(199) 


A/here 


A = k i k 0 . 
r rl  r2 


In  (198)  it  can  be  shown  that 


(200) 


Tr(nTf+nrTa)  = Tr(af+nrV  + (n-^Tf-cx(n^)-Co(n^)  • (201) 

Using  (194)  and  (201),  we  obtain  from  (198) 

\ ATrf 

cc  . Ar[T0(«TfanrTa)-Tr(«V„rTa)]  - Ar(-f- - j Tf 
A V f 

* if  L V„(tTf+"rTa)-2To(nTftnrTa>tTc("Tf+nrTa>-£r("Tf+nrTa)+Co<n'‘) 


■ 


1 


Using  { 194)  and  (202)  in  (193) 


VW'Wil-V'rV'WV  ■ (l-VtyiVrvV 

- Tr(<tTf+nrTa)]+W1+W2  (203) 


where 


W1  = VWaJ-y'rVWal-V'1^) 


AT 


rf 


and 


N -1 

Vf  -M-T--* 


A V1 

+ *7  n=0  iTo(nVrlrTa)-To^VnrTa)-Co(r1’^! 


A V 

w2  = R;  Jo  fo(nTf+nrTa,-Tc(nTf+nrV+  «r(nTf+nrTa)) 


(204) 


(205) 


We  shall  assume  that  the  clock  loop  is  in  the  steady  state.  Then 
making  use  of  (178)  and  ( 1 84. ) - ( 1 86 ) we  obtain  the  expected  values  of 
(204)  and  (205), 


ECW-,  ] = 


E[T0(-NTTf)]Tf 
N 
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p£(p  -1 ) . 1l  v p -p 

1 - p 


A„  V A n 
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2(T-^)  + p • ~ - P (1-p) 


/m  A \ m /m  A ' 

(^7 + Tf)  + "7  (AS 


n2T2 
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2T 


N -1 

VMt 

L„ 


7aT  * AT  ,, 

J cf  + rf 

If  If  I T 


(206) 
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Crs(nlrV*Tf+nrTa>  - V'VrVI'-V"'1'-1'-"/1 


•I  N -1  \ /4T  . IT,. 

^ ^ ^ V*c,<‘W.>)  • 


(213) 


For  the  system  to  be  stable. 


lim  (1-Ar)n=0. 

n-x» 


(214) 


Thus  we  must  set 


0 < A <2. 
r 


(215) 


In  the  steady-state,  the  mean  ranging  loop  timing  error  is  given  by 


fl  N -1  VATrf  ATrA 

ers(nIrTf+tTf+nrTa)  = ' l^/T  ' ~ +l)[~  + T^J  Tf+ecs(ATf+nrTa)  * 


(216) 


The  relationships  between  E[T0],  E[TC],  and  E[Tr]  for  E[To]=0,  Nr=2,  and 
Nc=l  are  shown  in  Figure  21.  This  figure  will  be  examined  in  detail  later 
(in  Section  V I - B)  when  applications  to  system  design  are  considered. 

To  calculate  the  variance  of  the  ranging-loop  timing  error,  we  assume 
that  all  the  random  processes  of  (203)  are  zero-mean.  This  is  permissible 
since  the  control  system  is  approximately  linear  in  the  sampled-data 
sense  over  the  range  of  operation  for  which  the  model  is  applicable. 

Using  ( 1 84 ) - ( 1 86 ) , (208)  and  the  following  relations 

Tc(nTf)  = To(nTf)-Te(nTf),  T.(nTf)  = fQ(nTf )-t£(nTf ) (217) 

Tc(nTf)  = TQ(nTf)  - Te(nTf)  (218) 

we  can  rewrite  (204)  and  (205)  as 
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Figure  21.  The  steady-state  delay-function  responses  of  the 
modified  synchronizer  for  the  case  Nr=2,  N =1 . 
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(cont.  next  page) 
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j2  Ye(nTf )+Te(nTf )nrTa  + T£(nTf)  Cr(nTf+nrTa)^  . 

(220) 


Due  to  the  complexity  of  the  expression,  we  shall  consider  two 
cases  only.  For  the  first  case  let  the  clock  loop  and  ranging  loop 
measurement  noise  variances  be  identical  and  assume  no  maneuver  noise, 
i . e . , 
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In  this  case,  all  the  random  variables,  u and  T , t , T are 
Thus  (203)  reduces  to 
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For  the  second  case  we  assume 


nr  = mT  = 0,  Ar  = Nr  = 1,  p > .99,  4 = Nx  = Ir  » 1. 


Because  the  largest  variance  occurs  at  t=NT,  this  is  the  most 
important  case  to  investigate.  In  this  case,  because  the  value  of 
is  close  to  unity, 
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and  (203)  reduces  to 
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Then  substituting  t = N into  (226), 


ejl  Tf+N  Tf+n  T ) % - \ [(N-k)2u(I  +k-l  )+N2u(k-l  )+k(2N  -k) 

r r T T i r d |^=  -j  T r T T 


f 

u( -N  +k-l ) J -x — 


' Wf>  ' V'rWf 


+ T (-»  T.)  t T (-N  T,)N  T.  + T (0)  * ( (n  T ) 
c t f e t f t f e r r a 


(227) 


and  substituting  NT=Ir  into  the  right-hand  side  of  (227)  and  taking  the 


't  1 r 

9 

expected  value  of  we  obtain 


lr  t4 

At  (WN,Wa>  * J C ( I r-k)4+I J+k2(2 Ir-k )2]  ^ (l-»2)»2 


* + °5r 


+ ^ El(k-l)[(Ir-k)2Te(IrTf) 


+ W 1! 


(228) 


or 
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<rarVNtTftnrTa)«  l [(Ir-k)4tI4rtkZ(2Ir-k)"]  jIO-b^o 

k 1 


T 


+3P-.  + 4P9I  T,  + 2P9I2T^+o? 

1 2 r f 3 r f {r 


+ kh  [(Ir-k)SkT(Ir-kt,)T2f+Ir‘’at(Ir-k+,)IrTfl 


(229) 

where  use  has  been  made  of  (191)  and  (192).  In  applying  Equation  (191), 
terms  are  neglected  if  i > Ir. 

Graphs  of  the  ranging  loop  error  variance,  Equation  (229)  are  shown 
in  Figures  22  and  23.  If  the  propagation  delay  is  assumed  to  be  .25  sec. 
then  Figure  22  applies  for  Tf  = .01  sec.  and  Figure  23  for  Tf  = .001  sec. 

From  the  two  cases  considered,  we  can  obtain  an  approximation 

f°r  aer(IrTf+NTTf+nrV  when  Nr  > 1 by  reP1acing  cr2^  with  a2r/Nr  in 
Equation  (229). 

For  the  case  of  a maneuvering  terminal  and  stationary  satellite, 
derivation  of  the  ranging  loop  timing  error  follows  essentially  the  same 
procedure  as  described  above  with  the  following  changes:  The  equation 
for  the  effects  of  a propagation  delay  change  CT,  (Equation  (194)) 
is  replaced  by 


Ct(«,I  «)  ■ T0(IrTf«ynrTa)  - T0(»Tf+nrTa) 


+ To<IrTfMWa->  - To<1Tf+nrTa-') 


(230) 


and  the  gain  factor  A0  (Equation  (209))  becomes 


A0  = 0 (231) 

i.e.,  no  open  loop  (predicted)  velocity  or  acceleration  corrections  are 
needed  since  any  change  in  range  for  the  case  of  the  stationary 
satellite  is  immediately  evident  at  the  terminal. 
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In  the  steady-state,  the  expected  ranging  loop  timing  error  ers  is 
still  given  by  Equation  (216).  However,  the  equation  for  calculating 
the  variance  of  the  ranging  loop  timing  error  becomes 

M'rV^WV  ■ 0-V'r(‘VnrTa) 


+ (1-Ar) 


m Tr 
t f 


n,2  P 

* V»Tf- NT)  -4-  I2 

T T 2N  r 
c 


m T , 


T.(  W'Wf >+te (IrV‘ VN,Tf > -r1 


m2  P 

* Tt(IrV'VN,Tf)  T7  T? 

C I 


t V V 

N n^n 
r n=0 


m T, 


ill  1 r- 

Te(nTf-NtTf)+fE(r,Tf-NTTf) 


• a m r\ 

+ T (n-N  ) -V-T; 

T 2 IT  T 

c _ 


A V1 

+ TT  a 
"r  n=0 


n2!2 

T.(nTf)+tt(nTf)n  T tT  (nTf)  -fi 


+ 5 (nT,+n  T ) 
r f r a 


(232) 


For  the  case 


1 = Nt  1r’  Ar  = Nr  = T*  nr  = mT  = 0 


the  variance  is 


(233) 


‘’er^rV^WV  ^ 3^1  + alr 


(234) 


S)(esCofT/„°2dSst0  E<’uati0".<2»)  “f  the  stationary  terminal  case. 

Of  Pj/ag  vs.  rp  were  given  previously  in  Section  C (see  Figure  8). 
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SECTION  V 

UP-LINK  TIMING  - EXTENDED  ANALYSIS 


A.  Introduction  and  General  Analysis 

This  section  presents  a more  complete  analysis  of  up-link  timing 
error  including  development  of  a model  applicable  to  botn  the  stationary- 

satellite  moving-terminal  case  and  the  stationary- terminal  moving- 
satellite  case.  An  optimum  form  of  augmented  open  loop  coupling  from 
the  clock  loop  to  the  ranging  loop  is  developed  (needed  only  in  the  case 
of  the  moving  satellite)  which  differs  from  that  utilized  in  Section 
IV-F  in  that  only  an  acceleration  term  (Tp)  is  used  rather  than  a com- 
bination of  acceleration  and  velocity  (T0).  This  is  possible  since  a 
constant  satellite  velocity  (zero-acceleration)  produces  a constant 
offset  in  receive  clock  and  up-link  timing  which  is  compensated  for  by 
the  coarse  ranging  system  during  initial  lock-up.  Thus  no  tracking  is 
required  to  compensate  for  the  presence  of  a velocity  term  and  only  the 
acceleration  term  is  needed.  Elimination  of  the  velocity  term  should 
reduce  the  noise  coupled  from  the  clock  loop  to  the  ranging  loop  and 
hence  result  in  reduced  error  variance  under  noise-limited  operating 
conditions.  The  analysis  of  this  section  also  retains  additional  cross 
correlation  terms  not  included  in  the  previous  analysis  and  conse- 
quently should  produce  more  accurate  results.  Typical  data  and  compari- 
sons with  data  from  the  analysis  of  Section  IV-F  are  given  in  Section  VI. 
Additional  notation  used  in  the  subsequent  analysis  is: 

R(t)  'v  distance  between  the  relay  and  a user  terminal 

D(t)  'v  a delay  function  defined  as  D(t)  = R(t)/c 
where  c is  the  speed  of  light 

d,j(t)  ^ down-link  delay 

dt(t)  'v  round-trip  (total)  delay  = x 

du(t)  'v  up-link  delay  referenced  to  the  receive  time  base: 
du(t)  6 dt(t)  " dd(t) 

t D(n)'v  time  of  occurrence  (receipt)  of  the  n-th  ranging  pulse 
rp  in  a burst  of  Nr  successive  ranging  pulses 

Tpr  ^ ranging  loop  processing  delay 

Tp  ^ frame  duration  = NpTf  where  Tf  is  the  subframe  duration,  and 
Np  = I 'v  Number  of  subframes  per  frame. 
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Also  we  have 


ec(t)  = To(t)_Tc(t)  (235) 

er(t)  = T0(t)-Tr(t)  (236) 

and 

erc(t)  - ■ (237) 

It  is  assumed  that  the  clock  frequencies  in  various  parts  of  the  system 
are  related  by 

(238) 


(239) 


and 


f =f  =f  =f  +Af 
r c o a o 


f =f 
s a 


where 


fr  'v  ranging  loop  code  generator  clocking  frequency 
fc  ^ clock  loop  code  generator  clocking  frequency 
fa  'v  assigned  code  generator  clocking  frequency 
fs  'v  effective  code  generator  clocking  frequency  associated  with 
the  network  clock  signal  radiated  by  the  satellite  (relay) 

A model  of  the  two-loop  terminal  synchronization  system  with  simple 
cross  strapping  applicable  to  the  stationary  satellite  case  is  given  in 
Figure  24  and  an  illustration  of  the  delay  function  vs.  time  for  a typical 
terminal  maneuver  is  given  in  Figure  25.  From  these  figures  we  obtain 


dd(t)  = D(t)  (240) 

dt(t)  = dd(t)+,,du(t),,-D(t)+D(  t-2D(t))+6Q(t-2D(t))[D(t)-D(t-2D(t))] 

(241) 

T0(t)  = D(t)  (242) 
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Figure  24.  Model  of  the  TDMA  synchronizer  for  the 
stationary  satellite  case. 


and 


VVTF>-  Tr(t„>  + VW  - VV 


Af 


t VWVVTF»  - WW  - -T  TF 


Af 


"c'VWVV'  - Tc«to-dt«to»  + -rM 


, N 
A . r 


+ R7  j,  [Ic<V"»  ‘ Wn))  - ‘n,<trp(n»] 


where 


‘o  ' [tcr  * MV)]  ' ‘o  + TF 


V = Vp'M  + V 


and 


* 2 TF  ' "7  rjl  [£F(trP<"»  - MV"11  - WV"))]- 

(247) 


A Nr 
Z v 


(243) 

(244) 

(245) 


6t(t)  = W + D( t+2D( t) ) + D( t+2D( t) ) [D( t+2D( t) )-D( t) ] . (246) 

It  follows  that 

MW  ■ MV-  MWMW1  * MVMV1 

Af 


V n=l  ' r'P  ' c'  rp'  " ’rn'  'rp' 
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Similarly,  from  Figures  26  and  27,  for  the  stationary  terminal  case 
there  results 


d (t)  = U(t-D(t) ) + D(t-U(t))  [D( t) -D( t-D( t) ) ] , 


(248) 


dt(t)  = 2dd(t), 


T0(t)  6 D(t-dd(t)), 


(249) 

(250) 


and 


Af 


T (t  +TC)  = T (t  ) + 2[D(t  +Tc-d  ,(t  +T_))  - D(t  -d  ,(t  ))]  - T_ 
rvo  F r o o F d'o  F " - o d ' o f F 


Af 


- [T  (t  +TF-d.(t  +TC))  - T (t  -d  (t  ))  + Tc] 
cx  o F t o F"  c o t o f FJ 


A„  Nr 


+ fT  ^ [Tc(trp(n))  ' Tr(trp(n))  " trn(trp(n))] 


(251) 


where 


t < [t  + 6. (t  )]  < t +Tr 
o cr  t cr  o F 


= t (N  ) + Tn 
r rp  r Dr 


(252) 

(253) 


and 


«t(t)  = 2[D(t+D(t) ) + D(t+D(t))[D(t+D(t))  - D(t)]]  . (254) 


Since 


T0(t)  = D(t-dd(t)) 


(255) 


it  follows  that 
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Figure  26.  Model  of  the  TOMA  synchronizer  for 
the  stationary  terminal  case. 
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DELAY 


MW  = MM  ~ M WMW)+  MVMM* 


Af 


A..  Nr 


+ 2 -M  "^c  ~ TT~  1 ^e„(t  (n))  - e (t  (n))  - t (t  (n))] 

f_  F N._  r rp  cv  rpv  rn  rp 


r n=l 


a (t  +Tr ) + a (t  ) 
o o F o o 


(256) 


where 

AQ(t)  = To(t)  - T0(t-dt(t))  . (257) 

Then  by  noting  that 

Ao(t)  “ A0(t+6t(t))  (258) 

A^(t-dt(t))=  AQ(t-dt(t))+  6t(t-dt(t))  = AQ(t)  (259) 


where  the  prime  indicates  an  equivalent  term  shifted  to  the  opposite  end 
of  the  delay  blocks  (d^(t)  = dci(t)  + du(t))  the  equivalent  model  of 
Figure  28  is  obtained. 

The  effect  of  the  term  A0(t+6^.(t))  in  the  preceding  model  (Figure  28) 
can  be  minimized  via  additional  open  loop  corrections.  Assume  that  such 
corrections  are  made  every  Tf/Na  seconds  and  that  ta  represents  a 
correction  instant.  Ideally,  the  open-loop  correction  at  t=ta+Tf/Na_  would 
equal  Oa(ta+tf/Na")  where 


°.<VW>  ■ -(Aot(ta+VNa 


5 - VVW  ' st<  WNa'>  ' r ■ (260) 

d 

Since  T0  is  not  known  exactly,  the  added  open  loop  corrections  can,  at 
best,  be  estimates,  0a,  where 


+ vvvv 


)]  - W6t(ta)] 


(261) 
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Figure  28.  Equivalent  model  of  the  TDMA  synchronizer 
for  the  stationary  terminal  case. 


The  addition  of  this  term  via  an  augmented  cross  strap  from  the  clock 
loop  to  the  ranging  loop  is  shown  in  the  augmented  model  of  Figure 
29. 

An  equivalent  model  for  the  TDMA  synchronizer  with  augmented  cross 
strap  is  given  in  Figure  30.  Although  appearing  somewhat  more  complex 
than  that  of  Figure  29  in  that  there  is  a certain  amount  of  redundancy 
in  the  arrangement  of  signal  inputs  and  function  blocks,  it  is  more 
amenable  to  the  analysis  to  follow  and  is  applicable  (by  setting  appropri- 
ate terms  to  zero)  to  the  stationary  satellite  case  as  well  as  to  the 
stationary  terminal  case.  From  Figure  30  we  obtain 
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Figure  29.  Augmented  cross  strap  model  of  the  TDMA  Synchronizer 
for  the  stationary  terminal  case. 
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Figure  30.  Equivalent  augmented  cross  strap  model  of  the 
TDMA  syncrhonizer. 


T (t  +TC)  = T (t  ) + T (t  +TC)  - T (t  ) 
rv  o F7  r o o'  o F7  o o 


•VVVVW)  - WMV’  --rT F 


- 1WWVVI  - Tc<Vdt(tol>  * -rV 


A Nr 

+ rr—  I [T  (t  (n))  - T (t  (p))  - t (t  (n))] 
N _i  c rp  r rp  7 rn  rp 


+ Q(VTF)  " Q(to) 


where 


t + 6.  (t  ) - Tc  < t < t + <5.  (t  ) 
cr  t cr  F o cr  t cr7 


^cr  Vp^r^  + T0r 


6t(t)  = 2^D(t+D(t:))  + D(t+D(t) ) [D(t+D(t))  - D(t)] 

TF  > dt  + Tf  + TDr 

Q(t)  = T0(t)  - T0(t-dt(t))  - Ta(t-dt(t) ) 


(262) 


(263) 


(264) 


(265) 


(266) 


(267) 


and  Ta(t)  is  defined  in  the  equivalent  model  of  the  augmented  cross 
strap  synchronizer,  Figure  30.  Similarly, 
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Tr(trp(n))  > Tr(t^T0(trp(n))  - T0(t0) 


if. 


+ To(trp(n)-dt(trp(n>»-  To(Vdt(to»  - rHJ"!'!,1 


Af 


-iTc(trp(n)-dt(trp(n)))-  Tc(to"dt(to))  + ~r  Ctrp(n)-to] 


+ Q(trp(n))  - Q(t0). 


(268) 


Using  this  expression  and  the  defining  equations  for  ec  and  er,  Equations 
(235)  and  (236),  it  can  be  shown  that 


. N 
A r 
_L  v 

Nr 


nl,  [TcVn))  ’ Tr(Vn))]  = V%(t0)^c(Vdt(t0))+Q(t0)] 


r [£c(trp(n))+ec(trp(n)-dt(trp(n))+Q(trp(n)^ 


(N-D  2Af 

Ar  ^to'trp^^  2 Tf^  T 

d 


(269) 


Incorporating  this  result  in  the  expression  for  Tr(t0+Tp),  Equation  (262), 
and  solving  for  E^(to+Tp)  gives 

Er<VV  ’ 0-Ar)[cr(t0)*ec(to-dt(toSH-Q(to)3 
- tclWdtlVV)  - «VV 

, N 

'i;  [ec(trp(n))tEc(trp(n)-<lt(trp(n,^Q<trp(n))+trn(trp(n>>] 


, (N  -1  ) I 24f„ 

MAr[(V‘n.(,»  --T-  V + TF 

a 


(270) 
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From  this  expression  it  is  evident  that 


er(to+2TF>  = ('-\)[V(totTF>«c(to+TF-dt(VTF>)t0(VrF» 

- ec<V2TF-dt<V2TF»-(’<V2TF> 


. N 
A r 
r v 


+ r n^/£c(trp(n^TF)+ec(trp(n)+TF-dt(trp(n)+TF^ 


+ Q(trp(n)+TF)+trn(trp(n)+TF^ 


+ <Art(VTF-(trp(l)+TF))  ' 


(N  -1) 


Tc]+T 


2Af 


FJ  F/  f 


(271) 


Substituting  the  expression  for  er(t0+Tp),  Equation  (270),  into  Equation 
(271 ) gives 

'r<V2V  ■ <'VK<V+ec(V‘VV^<V] 

- ec<V2Vdt(V2V>  - «V2V 


1 A ^r 

+ L (1-a/  / l [e_(t  (n)+(l-A)TF) 

4=0  r r n=l  c rp  h 


+Ec(trp(n)+{1-a)Vdt(trp(n)+{1-A)TF)) 

+Q(trp(n)  + (l-OTF)+trn(trp(n)  + (l-A)TF)] 


+<Ar[(t0-trp(l)) 


(N  -1) 


T,]+T 


2Af  1 


fJ  "V  f 


a ^=0 


( 1 -A  )‘ 
r' 


(272) 
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Continuing  in  this  manner  it  can  be  shown  that 


£r(to+mTF>  = 0-Ar)ni[er(to)+eC(to"dt(to))+Q{to)] 

- ec  ( to+mVdt  ( t0+mTF  ) ) "Q  ( t0+mTF ) 

m-1  A 

+ /V  ^"Ar^  if  ^ [e  (t  (n)+(m-l-4)T  ) 

4=0  r \ n=l  c rp  F 

+ ec^trp(n)+m"1"i)TF*dt(trp(n)  + ^'1"i^TF)) 

+ Q(trp(n)+(m-l-«.)TF)+trn(trp(r.)+(m-l-A)TF)] 


+ <V(Vtn,d))- 


(N  -1) 


T<J  + Tc 


2Af  m-1  „ 

l (1-A  )f 
Ta  4=0  r 


(273) 


Clearly,  (1-Ar)m  must  diminish  to  zero  as  m is  increased  to  a large  value 
if  the  transmit  timing  loop  is  to  respond  in  an  acceptable  manner. 
Consequently,  a necessary  condition  for  a "stable"  response  is  that  Ar 
assume  a value  between  zero  and  two: 

0 < Ar  < 2 . (274) 

Assume  that  this  condition  is  satisfied  and  that  m equals  a (any)  number 
M for  which 

(l-Ar)M«<l.  (275) 


It  follows  that 
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w 


and 


err,(to*TF>  S er(to‘*flTF3  " Ers(to*TF> 


* Jo  "-Ar>tfS<'cn'trp'M-1-1)  + S<‘c„.trp.M-l-‘) 


+ s<v„-vn-,-‘)] 


where 


£cn(t)  =ec(t)  ' E{£c(t)}  * 


(282) 


(283) 


It  follows  that 


= 4vS<V"TF>>  = ^ ^ W'V 


M-1  M-l  (i  +£  ) f 

+ jyo  if 0 1_Ar)  r^cn’Vp’^'V^cn’Vp’^'V 

* E{S<Ecn.VdfK-,-h)S<Ec„-Vp-dfM-'-12)} 

* E{S(trn’trP'H-1-tl)S(tr„-trp’M-,-‘2)| 

+ 2E{S<Ecn'trp’M-’-‘l>5<Ecn-Vdt'M-,-t2W 

mi  r /-  % — 


' 2 z=V1_V  Hcn(to^TF-dt)S<ecn*trp*M-1-A) 
+ E{cn<to*TF-dt)5<Ecn-VdfM-’-‘j>- 


(284) 


r 

i 


It  follows  from  the  defining  equation  for  S(x,t  ,k)  (Equation  (278)) 
that  rP 


ESS(tm*VM-1-V  S(trn’trp’M"1:i2) 


A s2  Nr  Nt  r 

[f)  lu}  E{Vnttrptnl)+<M-1-Jll)TF)trntVn2)  + (M-1-t2)TF) 
r 1 2 L > 


r 2 

N °t 
r rn 


*'2_il 


V*i 


(285) 


M-l  M-l  K+O  j 

I I (1-A  ) E/S(t  t ,M-l-4  )S(t  ,t  _,M-1-*  ) 

t1=0  t?=0  r rn  rp  I rn  rp  2 


Ar  2 Mr] 


rf  I ( 1 -A  )“  = 

r trn  4=0  r 


(286) 


r'  r rn 


Neither  ecn(t)  nor  ern(t)  are  stationary  in  the  strict  sense; 

consequently,  a?  (t0+MTp)  depends  on  (t0-dt),  t (1),  and  t (l)-dt)  in 

addition  to  other  parameters.  However,  the  expression  for  the  variance 
of  ern(t0+MTp)  is  reasonably  tractable  when 


trpO)  - (N1+ai)Tf 


(287) 


dt  = NdTf 


(288) 


tQ  = [No+(N1+a1)+Nr+N(J-NF]Tf 


(289) 


where 


0 < a.|  < 1 


N]  and  Nd  represent  constant  integers,  Nq  is  any  integer  for  which 


and 


0 I Nq  < Np  - 1 


0<To  <Tf 

r 1 


(290) 


(291) 


(292) 


Designate  the  variance  of  ern(t0+MTF)  when  these  constraints  apply  by 

%n  (rW  and  assume  that  clock  loop  has  been  operating  over 
an  indefinitely-long  interval  of  time  so  that 

E|cn(tl+KlTf)£cn(t2+l<2Tf^  = E(£cn  V VK)  VEcn(  V VK)  v} 


(293) 


where  t]  and  t 2 represent  arbitrary  instants  in  time  and  Ki , K?  and  K 
are  arbitrary  integers.  Finally,  define  (K.a)  where  I’is  an  ?„te 
as  ecn 


Rsc„(K’“>  ^ E{Ecn«K*°>Tf>scn(“Tf)  • 


nteger 


(294) 


Incorporating  the  notation  and  constraints  presented  in  the  Drecedino 
paragraph  into  the  expression  for  ^(t^Ty)  (Equation  (m))g?,es  9 
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L. 


/A  \2  M-l  M-l  (t,+0Nr  Nr 

+ 2(]f ) X I 0-Ar)  2 I I tRe  ((t1-A2)NF-(n1-n2),a1) 

Vr/  jyo  n-j  = 1 n2»l  ecn  1 


+ Re  ( (^i  ”*’2^F”^d” ^1  ”^2 ^ ,al  ^ 


A..  M-l 


-2  rf-  l (1-A  )z  l [R  (AN F+N  +N  -(n-l),ai) 
r a=o  n=l  cn 


+ Re  (ANF+Nd+Nc+Nr-(n-l),a1)].  (295) 
cn 


Using  this  equation,  it  can  be  shown  that 

A 


2A  Nr_1 

l (N  -ln|)[R  (M.O+R  (iN.+nj.a,)] 

(2-Ar)Nr2  |n=-(Nr-l)  r ecn  1 ecn  d 1 


+ I (l-Aj4  I (N  -|n|)[2R  (ANp+n.a.) 
A=1  r n=-(N  -1 ) r ecn  1 

"S 


+ R£  (|ANF-Nd+nl,a1)+Re  (ANp+Nj+n,^)] 


2A  « V 

ir  l l [Rc  (ANp+N  +n,oi  )+R£  (ANF+Nd+N  +n,ai)]. 

r a=0  n=l  cn  cn 


(296) 


When 


°1  <<  1 and 


2 

°e  (No»®l>  is  designated  as  a 2 (N)  and  is  given 


rn 


ern  0 


by 


\n(V^i+OT-‘t2 

rn  r'  r rn 


2A 


(2-A  )N  L 
r r 


V 

v 


(Nr-1 ) lNr' 1 " 1 1 CP1 ( |n I )+P1 ‘ 1 Md+n ' > ] 


V 


+ P^KNp-Nj+nD+P^tNp+Nj+n)] 


2A. 


Nr  Jo0"^^  nHf5l(ANF+Vn)+^^NF+Nd+No+n^*  (298) 


Similarly,  when 

0"°])  « 1 and  (1-^)  < 


(297) 


(299) 


°e  1S  desi9nated  as  a (N  ) and  is  given  by  Equation  (298) 

rn  rn 

replacing  ^(k)  by  ^(k). 


on 
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C.  Stationary-Terminal  Maneuverinq-Satel 1 i te  Case 

Now  consider  the  case  where  the  terminal  is  stationary  and  both 
normal  and  auxiliary  cross-straps  are  implemented.  Computational 
difficulties  are  minimized  by  presuming  that 

A =N  =N  =1 
r r a 

(300) 

Vp<'>*<V°i>Tf 

(301) 

VWf 

(302) 

and 

Vtv'W’V'-MTf 

(303) 

where 

0<  ( 1 -a-j  ) « 1 

(304) 

1 

P 

A 

(305) 

N]  and  Nj  are  constant  integers,  N0  is  any  integer  for  which 

0 1 N0  1 (Np-1) 

(306) 

and 

0 < Td  < Tf  . 
r 

(307) 

The  sampler  associated  with  the  auxiliary  cross  strap  is  presumed  to  be 
actuated  at  instants 

ta(K)  = KTf*TD+«aTf 

(308) 

where  K is  any  integer,  TD  represents  the  clock  loop  processing 
and  aaTf  is  a positive  lag  introduced  to  avoid  the  simultaneous 
of  two  series-connected  samplers.  It  is  presumed  that 

delay, 

closure 

TD  ~ Tf 

(309) 

and  that 

0 < a « 1 . 

a 

(310) 
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When  the  constraints  delineated  in  the  preceding  paragraph  apply, 
q^to,trp’M^  (E9uation  (28°))  is  designated  as  q^N^N-j.M)  and  is  given 

9!  ,M)  = QC(N0+Nd+l+N1+a1+(M-l ) Np ) Tf ]-Q[ ( N] + (M- 1 )Np)Tf]. 


by 


(311) 

Using  the  defining  equation  for  Q(t)  (Equation  (267)),  it  can  be  shown  that 
Q[(No+Nd+1+Nl+al+(M-1  )NF)Tf]  = Q[(N1+a1+(M-l)Nf:)Tf] 

+ To[(V1+VrV1+(M-1)NF)Tf>To[(Nl+1+(M-1)NF)Tf] 

" fo[(No+1+Nl+1+(M"1)NF)Tf]_Tol:(Nl+1"Nd+(M'1)NF)Tf]) 


Nq+1 +Nd 

°j1 


since 


T0[(K+ai)Tf]  = T0[(K+l)Tf] 


when 


0 < ( 1 — ct-j ) « 1 . 


(312) 


(313) 


(314) 


Using  the  foregoing  expressions,  it  can  be  shown  that 

To[(V1+Nd+Nl+1  + (M-1)NF)Tf]  = y(VUNl+U(M"1)NF)Tf] 

+ NdVoC(V1+Nl+1+(M_1)NF)Tf] 


+ 21  Sp(Nd)V(No+1+Nl+1+(M'1)NF)V 

JZf  Nd 

+ l Sp(Nd-k)u[N0+l+N1+l+k-l+(M-l)NF] 


where 


2K(1-p)-(1+p )(l-p  j 

A (1-P)2 

SP(K)  " 

0 , 


K > 0 


K = - 1 . 


(315) 


(216) 


The  latter  definition  makes  the  analysis  applicable  when  Nd  equals  zero. 
Similarly, 

T0[(N1+l+(M-l)NF)Tf]  = T0[(N1+l-Nd+(M-l)NF)Tf] 

+ NdTfto[(Nl+1-Nd+(M‘l)NF)Tf] 


+ r Sp(Nd)T0[(N1+l-Nd+(M-l  )Np)Tf] 


t!  Nd 


+ l Sp(Nd-k)u[N1+l-Nd+k-l+(M-l)NF] 


(317) 
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Combining  the  above  expressions  in  an  appropriate  manner  gives 


q1(N0,N1,M)  = NdTf<to[(No+l+N1+l+(M-l)NF)Tf]-to[(N1+l-Nd+(m-l)NF)Tf] 


+ 2~  Sp(Nd)jro[(No+l+N1+l+(M-l)NF)Tf]-To[(N1+l-Nd+(M-1)NF)Tf] 
T2  Nd  ^ 

+ J-  } S (N,-k)  <u[N  +l+N,+l+k-l+(M-l)NF] 

k-’  l ° , 

-u[N1+l-Nd+k-l+(M-l)N|:]> 


V1+Nd  r 

-Wt  J,  if0[(NrVkt(H-')NF>Tf] 


"Te[ta(NTNd+k+{M_1)NF)]/  (318) 


whenever 


T [t  (K)]  = T (KT  ) - T [t  (K)]  . 

c a o t e a 


(319) 


Tq[ (Nq+1  +N]  +!  +(M-1  )NF)Tf>t0[(N1+l-Nd+(M-l  )Np)Tf] 


(V1+Nd} 


Tq[ ( ni +T - Nd+( M-l  )NF)Tf] 


N +1+N.  f (N  +1+N  ,-k) 
0 V d l-p  0 d 


u[N-,  +1  -Nd+k-l  +(M-1  )Nf] 


(320) 


K +1+NJ  .. 


To[(No+1+Nl+1+(M"1)NF)Tf]  = p ° d 


V^d  (No+l+Nd-k) 

+ L P 

k=l 


ToC(Nl+1'Nd+(M'1)NF)Tf] 


u[N1+l-Nd+k-l+(M-l )Np]  . 


Since 


T0C(n+l )Tr]  = t0(nTf)+Tf.To(nTf) 


(321) 


(322) 


it  follows  that 


T0[(n+*)Tf]  = to(nTf)  + Tf  ^[(n+k^T^  , % >1. 


(323) 


It  was  previously  shown  that 


T0L(n+*)Tf]  = to(nTf)+Tf|^_^To(nTf)  + ^ u(n+k-l)|  , 


4 > 1 . (324) 


Jo  V(ntW  • (t^)  V"V  + j 0^-)  »<-M)  , 


4 I 1 . (325) 


Combining  the  last  few  equations  gives 


w 


1/  {Vund]\ 

ql(N0,NrM)  = NdTf\|~ To[(N1+l-Nd+(M-l)NF)Tf] 


V1+Nd  f,  <V1+Vk> 
+ y hz 


L 

k=l 


u[N1+l-Nd+k-1+(M 


Tf  (N  +1+N.) 

*rVV  >(’  d 


-1)  T0[(N1+NNd+(M-l)NF)Tf] 


V;Nd  (No+l+N  -k) 

+ £ p u[N1+l-Nd+k-l+(M-l)NF] 


T?  Nd 


f / 

+ T ^ Sp(Nd-k)\u[No+1+Nl+1+k-1+(M'1)NF^ 


- A N .H< 
a d f ' 


u[N1+l-Nd+k-1+(M-1 )Np] 
(No+l+Nd)' 


1-P 


TTP T0t(Ni+1-V(M-1)NF)Tf^ 


V1+Nd  L (V1+Nd"k^ 

+ y / 1-p  ' 


L 

k=l 


1-p 


/ 


u[N1+l-Nd+(M-l)NF; 


N0+1+Nd 

tAaNdT5  J,  Vta<NrNd+kt(M-1)NF>]  . 


(326) 
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1 


T2  Nd 

VwM)  :rU,  [sp(Hd>-sp<Vt>]u|;"rVlt+(H-1)NF5 


N +1+N  . 
o d 


+ S (N.)  l u[NrNd+k+(M-l)NF] 
p 0 k=l+Nd  1 d F 


+ I Sp(Nd-k)u[N0+l+N1+k+(M-l)NF] 


+ ( 1 +p  ) 


Nd\V1+Nd 


T£[ta(NrNd+k+(M-l)NF)]  . 


For  the  model  under  consideration, 
E {u(K) } = 0 


(329) 


(330) 


E {T£(KTf)}  = 0 


As  a result. 


wvNi’H)r°  • 


(331) 


(332) 


Designate  the  variance  of  q,  (N  ,N,  ,M)  by  aj-  (N  ).  Using  the 
expressions  p 0 qlp  0 


°u  ’ K1  = K2 


E/u(K1)u(K2)>  = 


Eye[ta(KlTf)]  u(K2J  = 0 • K2  > Kl 

and  the  foregoing  equation  for  q1  (Nq ,N1 ,M) , it  can  be  shown  that 


(334) 
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”L<No>  s 


_H(Sc(Nd)-sp(Nd-k))V(N(J-k)]+(No+l)Sp2(Nd)LE 


NdA  N0  + 1+Nd  V1+Nc 


y 

Z 


v 

L 


k!=l  k2=l 


f Te[ta(NrNd+kl+(M-^NF^ 


‘ Te[ta{NrNdV(H-,)NF)] 


f Nd\l  N0+1+Nd  kl 

+ 2(,+B\«C<N.)  J.,  kI=,  EjTe[ta(NrVk]+(M.1)„F)] 


T 


* utNrNd+k2+(M-l)NF] 


V+Nd  Nd 


J=1  k 1=1  SP(Nd’k2)E^[ta(NrNd+k1+(M-1)NF)] 


‘ u^rNd+k2+(M-l)NFA>}  . (335) 


Equivalently, 


V 


°qlp(No)  ^ TT  \fo*™**>J-*  JQ  SP^W~SP^V  °u 


+ (1+p) 


2 1-p 


T~p  I ^<V,+Nd>Ri  ‘“'“a  + T7> 


+ 2 ( 1 +p  ) 


No+Nd  i 

* 2 j,  (VV^C'V  17>! 


/ MVd 

1-p 


V 


its-  z rtu (*.«a+  f) 

P / 4=0  e a ‘f  \ 


where 


• tSp(Nd)(No+l+Nd-4)-Fs(4)] 

R-j  ( K,a)  = E^E[(K+K1+a)Tf]  Te [ ( ^ +a)Tf ]j 
Rj  u(K,a)  = E^£[(K+K.|+a)Tf]  u(K^)|  , 0<a<l 

and  are  arbitrary  integers,  and 


t n 


V 


l . WM  ■ I S» 


k0=l 


FS(M  * 


Nd-[4-(No+l] 
l S (Nd-k2) 


k2=l 


k=0 


V 

l S (k) 
k=4-(N%l) 


» 4 - V1 


i > N +1 
— o 


(336) 

(337) 

(338) 


(339) 
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The  random  component  of  the  ranging  loop  timing  j^rror  for  the 

particular  case  under  consideration  is  designated  as  e (t).  From 

rn-, 


earlier  expressions,  it  can  be  shown  that 

e'rn,  <VV'  - 


'Ip 


IP 


* ec„[(N,+a,t(M-l)NF)Tf] 


- VW"> 


provided 

0 < (1-a-j ) « 1. 

It  follows  that 

(N  ) = a2  (N  ) 

rnlp  ™ 


VV 


+ % <N0>  - 2\  * <No> 
Mlp  Mlp  rn 


where 


Vrn<V°  ' ' 4.P("o-VH)Ec"[(V'*W<»-’I»F>V 

- e cn [ ( N i +<x i + ( M- 1 )Np)Tf] 


• ecn[(Nl+arNd+(M-1)NF)Tf] 


(340) 

(341) 

(342) 


(343) 


Using  the  foregoing  equation  for  (N  ,N,  ,M) , Equation  (329)),  it  can' 
be  shown  that  ip  o i 
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HI 


N +1+N. 
o d 


SP(Nd)  kl}  Tcn^V’+VV^VV 


v 


■u[NrNd+k+(M-l)NF] 


- J1Sp(Nd-k)E|ecnC<No+1+Nl4al+^M-1iNF)Tf] 

•u[NrNd+k+(M-l)NF]' 

Nd  r 

• W J^ecn[(Nl+V(M-1)NF)Tf] 


■u[NrNd+k+(M-l)NF] 


•u[NrNd+k+(M-1)NF]j 

f Nd  \ Nq+1 +Nd  ( 

+(1+p)te“ ) kl}  EyE[ta(NrVk+(M-1)NF)] 

•[cnC(V1+NlMl*(H-,)NF)TfJ 

-tcn[(N,t»,*(H-l)NF)Tf] 

-tcn[<Nit“rV<M-,>NF>Tfi]! 


(344) 
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Equivalently, 


N +N  . 
o d 


-2R*  * (N  )n  = lUs  (N.)  I R-v  (t) 
^lp^rn  00  f 1 p d 'Zu  e™u 


*.=N,  cn 


V 

+ Jo  SpU)  \nu(M  ' \nu(UV1)] 


+ (1+p) 


M -p 


1-P 


N +1 
o 


'N  +N  . 
o d 


i t y 

£=N,  Ecn  e 
d 


N +1+N  , 
o d 


- 2 l Rj  y (*)  - I K y (*) 
*■=1  e£cn  i=N  +2  e*cn 


where 


R^cnu(K)  = E\ccn[{l<+l<l+al)Tf]  u[Kl]' 


- E(ecn[(K+Kl+al)Tf]  TeC(K1+a2)T  f]j 

RbjK>  « ecn[(K,t«])d 


e cn 


0 < ( 1 -a-j ) « 1 


0 < a2  « 1 


and  K and  K-j  are  arbitrary  integers.  Note  that 
R^  :j  HO  = R^  ^ (K). 


cn  e 


e cn 


(345) 


(346) 

(347) 

(348) 

(349) 

(350) 

(351) 
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TJe  equations  developed  in  this  section  are  applied  in  Section 
h?  SaiCJ!late+tyPlCal  data  aPPlicab1e  to  up-link  synchronization 
system  bit-synchronous  TDMA  satellite  (relay)  communications 


no 


SECTION  VI 

APPLICATIONS  TO  SYSTEM  DESIGN 


\ 


In  this  section  the  steady  state  solutions  for  the  clock  and  ranging 
loop  timing  error  variances  derived  in  Sections  IV  and  V are  applied  to 
determine  system  specifications  and  performance  as  a function  of  known, 
or  chosen,  system  parameters.  Since  the  form  of  the  solutions  obtained 
earlier  is  not  particularly  convenient  for  this  use,  different  forms 
are  derived,  and  design  curves  are  given  which  permit  optimum  system 
design  with  minimum  effort. 

In  Section  V I -A  only  the  clock  loop  timing  required  to  receive 
properly  timed  signals  on  the  down-link  is  considered,  while  in  Section 
VI-B  the  ranging  loop  timing  error  limits  adequate  for  transmission  are 
considered.  Because  of  the  coupling  betv/een  the  clock  and  ranging  loops, 
as  discussed  earlier,  an  adequate  design  for  transmission  timing  normally 
results  in  more  than  adequate  receive  clock  loop  performance.  Hence  for 
terminals  involving  both  transmission  and  reception  the  results  of 
Section  VI-B  are  applicable,  while  the  results  of  Section  VI-A  yield  con- 
siderably less  critical  specifications  for  receive  only  terminals. 


Clock  Loop  Timinc 


A typical  design  problem  would  be  to  determine  the  system  parameters 
required  to  establish  communications  at  a given  chip  rate  and  specified 
maximum  error  rate  with  a terminal  whose  maneuver  parameters  are  known 
or  can  be  estimated.  A method  for  system  design  utilizing  the  steady 
state  results  of  Section  IV-C  for  a receive-only  terminal  employing  an 
augmented  Kalman  filter  will  now  be  illustrated  by  example. 


Example  1 


Given: 


Chip  Rate  (i)  = 40  Mbit/s 

-4 

Maximum  Bit  Error  Probability  = 10 

2 

Maximum  Terminal  Acceleration  (A)  = 49  m/sec  (5g.) 
Probability  That  Terminal  will  not  Accelerate  (P  ) = .5 
Average  Maneuver  Duration  (1/A)  = 5 sec. 


Ill 


Find: 


Suitable  System  Parameters: 

Tf  = KMA  = subframe  time  (i.e.,  time  between  system  clock 
samples) 

Ma  = clock  pulse  length 

E/Nq  = clock  pulse  energy  to  noise  density  ratio 
P/Nq  = received  power  to  noise  density  ratio 

The  Kalman  Gain  Constants  (Steady  State): 

Kj  = range  gain 
K-f  = velocity  gain 
Kj  = acceleration  gain 

Other  Related  Parameters: 

a = receiver  (measurement)  noise  standard  deviation 
a.  = maneuver  noise  standard  deviation 

a 

p = correlation  coefficient 

P-|  = filtered  estimate  (minimum)  timing  error  variance 
= predicted  estimate  (maximum)  timing  error  variance 

From  a previous  study  [1]  the  maximum  timing  jitter  resulting  in  no 
significant  increase  in  error  rate  is  approximately  .05a.  Hence  we 
choose  the  predicted  estimate  (worst  case)  error  variance  as 

= (.05A)2  = 1.5625  • 10'18  (sec2).  (352) 

From  comparisons  of  typical  values  of  the  error  variances  and 
P-]  vs.  the  measurement  noise  variance  a^2  in  Section  III,  we  choose  as 
a reasonable  value  for  o 


= ,05a  =1.25  • 10"9  sec.  (353) 

Then : 


Using  the  maneuver  model  of  Figure  6 and  Equation  (80),  with  the 
probability  of  maximum  acceleration,  Pm,  set  equal  to  zero*  yields  the 
maneuver  standard  deviation 

— .1/2 

/\ 

(1  - PJ  = — = — = 6.67  • 10'8  sec/sec2  (355) 

°_|  /6  c/6 

where  c is  the  velocity  of  light. 


'v  2 

Now  from  the  functional  relationship  relating  P-|/a?  vs.  rp  for  various 
values  of  p,  Figure  7,  and  the  defining  relations  (from  Section  IV) 


r 

p 


= r 

r 1+p ' 


and 


(356) 


»-l-ATf  (Tfll) 


(357) 


The  value  of  the  subframe  time  Tf  can  be  determined  (by  iteration)  in 
the  following  manner-  From  the  intersections  of  the  curves  of  Figure  7 
with  the  coordinate  Pi/o  2 = i.o  (from  Equation  (354))  the  values  for 
columns  (1)  and  (2)  of  Table  1 are  obtained. 

Table  1 


System 

Parameters  for 

the  Solution 

of  Example  1 

0) 

(2) 

(3) 

(4) 

(5) 

p 

rp 

Tf  (sec) 

r 

o,  (sec/sec  ) 

a 

.99 

0.94 

.05 

13.26 

1.51  • 

10"7 

.9 

6.8 

.5 

29.64 

6.75  • 

IQ'10 

.5 

14.5 

2.5 

25.11 

3.19  • 

IQ'11 

0 

16.2 

5.0 

16.2 

1.23  • 

IQ'11 

*$ince  the  relative  acceleration  between  vehicle  and  satellite  is 
dependent  upon  direction  as  well  as  magnitude,  there  is  no  reason  to 
believe  the  maximum  values  are  more  probable  than  other  values. 
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Corresponding  values  for  column  (3)  and  for  columns  (4)  and  (5)  are  then 
obtained  via  Equations  (357)  and  (356).  respectively.  Then  from  Table  1 
and  the  value  oa  = 6.67  • 10"8  sec/sec^  obtained  in  Equation  (355)  we 
see  that  the  solution  must  be  somewhere  in  the  range  .05  < Tf  < .5  and 
.99  > p > .9.  Accurate  interpolation  directly  from  Figure  7,  however, 
would  be  very  difficult  because  of  the  few  (4)  applicable  data  points 
available  and  the  non-linear  spacinq  of  the  curves.  To  obtain  additional 
accuracy  it  is  noted  that  the  data  of  columns  (1)  and  (2)  of  Table  1 
yield  a smooth  curve  which  can  be  readily  interpolated  if  plotted  in 
the  form  rp  vs.  (1-p)  on  a log-log  plot  as  shown  in  Figure  31.  Additional 
values  for  Table  1 are  then  computed  for  values  of  p and  rp  obtained 
from  this  curve,  converging  on  the  required  value  of  oa  until  the 
desired  accuracy  is  obtained.  The  result  is  found  to  be 


r = 1.29 

P 

T^  = .0695  sec 


r = 15.42 

o = 6.7  • 10-8  sec/sec2  = 2.68a  (sec-2) 


(358) 


Then  from  Equation  (77): 


1.215 


sfe  Wir) 


1.215 


1.476C 


4( a ) 


(359) 


where 


2 for  square  law  detection 
1 for  linear  detection 


(360) 


and  substituting  the  values  obtained  above  (with  C=l) 


= 122.7  * 20.89  dB. 


(361) 


If  we  now  choose  the  number  of  slots  per  subframe 


K = 32  (as  in  the  present  system) 


(362) 


Then  the  timing  pulse  (slot)  length  becomes 


MA  = Tf/K  = 2.172  • 10"3  sec 

and  the  required  received  power  to  noise  density  ratio  is 
E/N  . 

P/Nq  = -j~  = 5.649  • 10  Hz  'v  47.5  dB  Hz  . 


(363) 


(364) 


The  Kalman  gain  factors  (steady  state)  needed  to  implement  the 
filter,  Ky,  Kf  and  K’j,  and  the  standard  deviation  of  the  timing  error 
just  after  a data  point  is  received,  (Py )3/^,  can  be  found  from 
Figures  8,  9 and  10  for  the  values  of  p and  rp  as  determined  above. 
Difficulties  in  interpolation  are  encountered  here  also,  however;  and  so 
auxill iary  curves  of  Kt,  KfTy  and  KyTf^  vs.  (1-p)  are  first  drawn  using 
points  obtained  from  the  intersections  of  the  given  curves  with  the 
coordinate  rp  = 1.29  as  obtained  above.  The  resultant  curves  are 
given  in  Figure  32  (together  with  similar  curves  for  example  2,  below). 
From  these  curves,  for  (1-p)  = .0139,  the  values 


Ky  = — y = .497 

°C 

K|Tf  = .167 

K"Tf2  = .029  (365) 


are  obtained  from  which 


Kt  = .497,  Kj.  = 2.40  (sec-1),  = 6.00  (sec-2) 


and 


(P^172  = 8.812  • 10"10  sec  = .035a  . 


(366) 


(367) 


Example  2 

This  example  is  the  same  as  Example  1 except  the  Bit  Rate  has 
been  increased  to  80  Mbps.  Results  are  obtained  as  for  Example  1 
utilizing  the  auxilliary  curves  of  Figures  31  and  32.  Parameters  for 
this  example  are  compared  to  those  of  Example  1 in  Table  2.  It  is 
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( I — /O  ) 

Figure  32.  Kalman  gain  factors  vs.  (1-p)  for  examples  1 and  2. 
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interesting  to  note  that  while  a decrease  in  subframe  time  (i.e.,  the 
time  between  network  clock  pulses)  and  an  increase  in  received  power 
to  noise  density  ratio  are  required,  as  expected,  for  the  higher  code 
rate,  the  relative  changes  in  these  parameters  are  rather  small 
compared  to  the  relative  change  in  code  rate.  Hence  it  appears  that 
maintaining  adequate  receive  clock  timing  accuracy  while  increasing 
the  code  rate  within  reasonable  limits  should  not  be  too  difficult  when 
a Kalman  filter  is  employed  in  the  clock  loop.  The  large  number  of 
code  chips  per  slot  for  both  examples  permits  long  codes  to  be  used 
effectively  and  a high  degree  of  spectrum  spreading  to  be  employed  if 
desired,  while  maintaining  reasonable  signaling  rates. 


Table  2 

SYSTEM  PARAMETERS  FOR  CLOCK  LOOP  SYNCHRONIZATION 


Parameter 

Examples  1 & -3 

Example  2 

Code  Rate  (1/a) 

40  Mbps 

80  Mbps 

Max.  Timing  Error  Std.  Dev.  (f^i)1//2 

.05A 

(same) 

Measurement  Noise  Stn.  Dev.  (of) 

.05a 

(same) 

Max.  Terminal  Acceleration  (A) 

5g. 

(same) 

Maneuver  Rate  (X) 

0.2  (Maneuvers/sec) 

(same) 

Maneuver  Noise  Std.  Dev.  (aa) 

2.68a  (sec'2) 

(same) 

Correlation  Coefficient  (p) 

.9861 

.9896 

1 o 1/2 

1.29 

.997 

4o 

r = -XT 

15.42 

13.789 

Vf 

Subframe  Time  (Tf) 

.0695  (sec) 

.052  (sec) 

E/N0 

20.89  dB 

(same) 

Slots/Subframe  (K) 

32 

(same) 

Timing  Pulse  Length  (Ma) 

2.1 72 • 1 0 ~ 3 (sec) 

1 .625- 10~3  (s 

Code  Chips/Slot  (M) 

86,880 

130,000 

P/N° 

47.5  dB  (Hz) 

48.78  dB  (Hz) 

Min.  Timing  Error  Std.  Dev.  (P-|)^2 

.035a 

.0351A 

Kj(range  gain) 

.497 

.493 

Kj  (velocity  gain) 

2.40  (sec-1 ) 

3.26  (sec-1) 

Kj  (acceleration  gain) 

6.00  (sec-2) 

10.17  (sec-2) 
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An  iterative  procedure  is  required  in  solving  the  above  examples 
because  the  principal  unknown,  Tf,  appears  implicitly  in  both  the 
parameter  rp  (see  Equation  (356)),  the  abscissa  of  Figure  7,  and  in 
p(see  Equation  (357))  which  identifies  a particular  curve  of  the 
family.  This  difficulty  can  be  avoided  by  calculating  and  plotting 
the  Kalman  gains  and  error  variances  as  a function  of  a new  variable 
r$  which  contains  only  known  parameters 


(368) 


Also  it  is  convenient  to  specify  each  curve  of  the  family  by  the 
parameter 

ATf  = 1-p  (from  Equation  (357)) 


(369) 


since  this  leads  to  more  uniform  curve  spacing  and  facilitates 
interpolation  as  noted  earlier. 

Since  the  original  solution  (Equation  (118))  was  obtained  as  a 
function  of  the  variables  p and  rp , to  obtain  the  desired  curves  an 
expression  relating  the  new  variable  rs  to  rp  and  p is  needed.  This  is 
obtained  by  combining  Equations  (356),  (368)  and  (369)  resulting  in 


r 

s 


(1-c)2 

4 


r 


(1-p)2 

4 


(370) 


Curves  of  the  Kalman  gains  and  error  variances  vs.  rs  are  given 
in  Figures  33  through  36.  Example  3 below  illustrates  use  of  these 
curves  in  solving  the  system  design  problem  of  Example  1. 


Example  3 

From  Example  1: 

(371) 

(372) 

(373) 
and 


o = .05a  =1.25  • 10"9  sec 
a = 6.67  ■ 10~8  sec/sec2 

d 

A = .2  maneuvers/sec 
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Figure  34.  Kalman  gain  and  filtered  estimate  error  variance  vs. 


Figure  35.  Normalized  Kalman  velocity  gain  vs. 


Figure  36.  Normalized  Kalman  acceleration  ga 


(374) 


= 1.00  . 

Then  from  Equation  (368): 

o A2 

r = = 7.50  * 10'4  . (375) 

a 

Using  the  values  from  Equations  (374)  and  (375)  as  coordinates, 
we  find  from  Figure  33  (by  interpolation  between  the  given  curves) 

ATf  % .014  (376) 

and  from  Figures  34-36,  (for  r$  = 7.50  • 10”4  and  ATf  = .014) 


-L  % .50 

(377) 

K|Tf  % .17 

(378) 

K-^Tf2  * .028 

(379) 

which  are  in  good  agreement  with  the  values  obtained  in  Example  1. 

Since,  as  noted  in  Example  1,  the  system  subframe  time,  Tf,  is  the 
principal  unknown  in  a typical  system  design  problem,  the  effects 
of  parameter  variations  upon  system  design  can  best  be  illustrated  by 
studying  the  effects  of  parameter  variations  upon  the  subframe  time. 
Curves  of  the  system  subframe  time  versus  maximum  terminal  acceleration 
are  given  in  Figures  37-41  for  a variety  of  system  parameters.  In  each 
figure  a specific  parameter  is  varied  (while  others  are  held  constant) 
to  obtain  a family  of  curves,  thus  showing  the  effects  of  this 
parameter  upon  the  resulting  subframe  time.  These  curves  were  obtained 
by  a computer  iteration  technique  converging  to  specified  values  of 
Pl/a?2  essentially  by  the  method  of  Example  3,  above.  The  values  of 
Tf  thus  obtained  represent  maximum  values  which  may  be  used  to  achieve 
the  specified  timing  jitter  within  limits  set  by  the  other  (fixed) 
pa rameters. 
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Figure  37.  Subframe  time  vs.  maximum  terminal  acceleration  for  different  values 
of  clock  loop  timing  jitter. 
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measurement  noise  standard  deviations,  clock  loop. 


From  the  data  of  Figures  37-41,  in  the  range  of  parameters  con- 
sidered, the  subframe  time  is  given  approximately  by 

Tf  % k f-j  A-*19  (1-P0P21  A-'41  A’42  (380) 

where  k is  a constant,  f]  is  a rather  complex  function  of  and 
and  the  other  symbols  are  as  previously  defined.  From  the  functional 
dependence  of  Tf  vs.  A (also  Table  2)  we  note  that  while  increasing  the 
code  rate  (decreasing  a)  requires  a decrease  in  the  subframe  time  as 
expected,  the  number  of  code  chips  per  subframe  actually  increases. 

This  increase  is  not  "free"  however,  since  if  the  extra  chips  are  to 
be  used  to  carry  information  (rather  than  simply  to  increase  the 
spectrum  spreading  ratio)  then  the  received  power  will  have  to  be 
increased  to  maintain  the  necessary  signal  to  noise  ratio  prior  to 
detection. 


B.  Up-Link  Timinc 


An  expression  for  the  mean  steady  state  up-link  timing  error  for 
the  stationary  terminal  (i.e.,  worst  case)  was  developed  previously  in 
Section  IV-F,  Equation  (216),  and  typical  examples  of  the  associated 
delay  function  responses  were  given  in  Figure  21.  This  figure  represents 
conditions  for  which  the  satellite  velocity  relative  to  the  terminal  is 
a constant,  r,  (range  increasing  with  time)  while  the  acceleration  is 
zero  (i.e.,  the  terminal  is  not  maneuvering).  The  rate  of  change  of 
propagation  path  delay  is  then  t0  = r/c  where  c is  the  velocity  of 
propagation.  The  subscript  s in  Figure  21  indicates  expected  values, 
i.e.,  the  effects  of  noise  are  not  included. 


Physically,  the  curves  of  Figure  21  can  be  interpreted  in  the 
following  manner.  The  curve  Tos(t)  represents  the  actual  path  delay 
versus  time,  while  the  curve  Tcs(t)  represents  the  clock  loop  (with 
Kalman  filter)  estimate  of  this  delay.  Hence  the  clock  loop  timing 
error,  less  the  effects  due  to  noise,  is  given  by  the  difference  in 

delay  between  these  two  curves  at  any  point  in  time.  The  mean  error  is 

seen  to  be  zero  due  to  the  prediction  capabilities  of  the  filter  and 

the  method  of  error  correction  used.  The  maximum  absolute  error  due  to 

the  discrete  corrections  in  clock  frequency  is  given  by 


and  hence  can  be  kept  small  (for  a given  clock  frequency  offset)  by 
Choosing  the  number  of  clock  frequency  corrections  per  sampling 
interval,  Nc,  sufficiently  large. 
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By  analogy  to  the  above  we  would  like  to  consider  the  curve 
Trs(t)  of  Figure  21  as  the  ranging  looo  estimate  of  the  propagation 
delay  but,  unfortunately,  this  is  not  the  case  since  an  actual  delay 
estimate  is  not  derived  in  the  current  ranging  loop  implementation. 

The  curve  Trs(t),  which  represents  the  timing  function  (sampling  due 
to  discrete  ranging  pulses  neglected)  at  the  point  Tr(t)  in  Figure  19, 
is  simply  related  to  such  a delay  estimate  however,  and  will  serve  as 
well.  By  considering  the  model  of  Figure  19  relative  to  the  physical 
situation,  it  can  be  seen  that  a ranging  loop  estimate  of  the  propa- 
gation time  delay  could  be  obtained  as 

T ,(t)  = T (t)  - (T  (t)  - T (t))  = 2 T (t)  - T (t)  (382) 

pd  os  rs  os  os  rs 

i.e.,  each  point  on  the  curve  Tpd(t)  is  the  vertical  reflection  of  the 
corresponding  point  on  Trs(t)  (same  value  of  t for  both  points)  about 
the  curve  T0s(t).  Hence  the  ranging  loop  timing  error  is  (except  for 
a sign  reversal)  given  by  the  difference  in  time  delay  for  the  curves 
Trs(t)  and  Tos(t)  at  any  time  t,  analogous  to  that  of  the  clock  loop. 

As  indicated  in  Figure  21,  measurements  made  of  the  ranging  loop 
timing  error  (relative  to  the  clock  loop  estimate  Tcs(t)  since  the 
actual  value  T0s(t)  is  not  available)  at  C and  D are  averaged  and 

a correction  made  at  D which  would,  in  the  absence  of  propagation 

delay,  superimpose  the  curves  Trs(t)  and  Tcs(t)  at  that  point  (within 
limits 'set  by  the  average  value  re  a five  to  the  value  at  D).  However, 
because  of  the  round  trip  propagation  delay,  x,  this  correction  is  not 
seen  in  the  curve  Trs(t)  until  x seconds  later,  at  which  time,  due  to 
clock  oscillator  offsets,  the  correction  will  not,  in  general,  be 
the  required  amount  to  superimpose  the  delay  functions.  Thus  the  mean 

ranging  loop  timing  error  is,  in  general,  non-zero  because  of  (1)  the 

effects  of  discrete  clock  loop  corrections  and  correspondi ng  open  loop 
corrections  to  the  ranging  loop  upon  ranging  loop  error  measurements 
(i.e.,  the  small  steps  in  Tcs(t)  and  Trs(t)  of  Figure  21)  and  (2)  the 
effects  of  propagation  delay  and  clock  oscillator  offsets  upon  the 
subsequent  closed  loop  error  corrections  (the  large  saw  tooth  waveform 
in  TrS ( t) ) . 

Note  that  the  open  loop  corrections  to  the  ranging  loop  are  delayed 
in  the  response  Trs(t)  relative  to  the  clock  loop  by  x seconds,  as  is 
the  closed  loop  correction.  Thus  the  alignment  of  the  offsets  in 
Trs(t)  relative  to  those  of  Tcs(t)  is  x dependent  and  could  take  on 
any  value  as  the  terminal -satell i te  range  is  varied.  If  we  assume  the 
error  measurements  are  made  near  the  times  nTf  (n  an  integer)  in 
Figure  21  so  that  the  error  between  Tcs(t)  and  Tos(t)  may  be  neglected 
(in  the  present  system  error  measurements  are  made  for 
(n  + l/16)Tf<t<(n  + l/4)Tf,  depending  upon  the  link  range  slots  assigned 
to  that  particular  terminal),  then  the  maximum  error  due  to  the 
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discrete  open  loop  corrections  (for  a non-maneuvering  terminal)  would 
be  given  by  Equation  (381).  For  non-optimum  sampling  times  the  maximum 
error  could  be  twice  this  amount. 


From  Figure  21,  the  maximum  ranging  loop  timing  error  due  to  the 
effects  of  propagation  delay  and  clock  oscillator  offsets  on  the 
closed  loop  corrections  is  found  to  be 


t.  rM 


»rTf 


(383) 


and  the  mean  error  due  to  this  cause  is 


Thus  to  reduce  this  error  (since  we  normally  have  no  control  over  the 
propagation  delay  t),  we  must  minimize  the  clock  oscillator  offsets 
(i.e.,  use  adequately  stable  clocks  in  both  satellite  and  terminals)  and 
keep  the  frame  length  IrTf  short,  consistent  with  other  requirements. 

Expressions  for  the  ranging  loop  timing  error  variance,  including 
the  effects  of  measurement  noise  and  terminal  maneuver  statistics,  were 
developed  previously:  Equations  (229)  and  (234)  of  Section  IV  and 
Equations  (298)  and  (342)  of  Section  V.  Assuming  the  mean  error  can  be 
kept  sufficiently  small  by  appropriate  choice  of  parameters  as  discussed 
above,  the  error  variance  becomes  the  crucial  design  criteria  as  for 
the  clock  loop  timing  of  Section  A (the  mean  error  for  the  clock  loop 
is  zero,  as  noted  above),  hence  application  of  these  variance  equations 
to  system  design  will  now  be  investigated.  The  development  essentially 
parallels  that  for  the  clock  loop  (down-link)  timing  of  section  A, 
except  that  an  additional  parameter,  Ir,  the  number  of  subframes  per 
frame  (which  together  with  the  subframe  length  Tf  determines  the  frame 
length,  1 rT f ) must  be  considered.  Also,  Nr,  the  number  of  ranging 
pulses  per  burst  to  be  averaged  in  determining  each  closed  loop 
correction,  is  assumed  equal  to  one  in  the  basic  analysis  with  cor- 
rections for  other  values  of  Nr  to  be  made  by  replacing  o?r2  with 
°£r  /Nr  as  noted  previously  in  Section  IV-F. 

Curves  of  normalized  ranging  loop  timing  jitter  versus  r for  the 
stationary  terminal  case  (from  Equation  (229))  were  given  in  Figures 
22  and  23  of  Section  IV.  However,  direct  application  of  these  curves  to 
system  design  is  not  convenient  because  of  the  implicit  dependence 
upon  Tf  of  both  r and  p as  noted  previously  for  the  case  of  clock  loop 
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timing,  Section  A.  Hence  curves  of  normalized  ranging  loop  timing 
jitter  vs.  rs  (see  Equations  (368)-(370)) , which  are  more  directly 
applicable,  are  given  in  Figures  42  through  44  for  different  values  of 
Ir.  Application  of  these  curves  to  system  design  will  be  illustrated 
by  an  example. 

Example  4 

Given:  Bit  Rate  (~)  = 40  Mbps 

Maximum  up-link  timing  jitter  (a  ) = .05a 

Measurement  noise  standard  deviation  (o ^ - a ) - .0375a 

Ranging  pulses  per  burst  (N^)  = 1 

Probability  that  terminal  will  not  accelerate  (P  ) = .5 

j ' o 

Maneuver  rate  (A)  = .1  maneuvers/sec 

2 

Maximum  terminal  acceleration  (A)  = 9.8  m/sec  (1  g.) 
Number  of  subframes/frame  (1^)  = 32  (or  greater). 


Then  t,.e  normalized  timing  jitter  is: 

/ 2 2X,  2 (.0 5a)2  - (.0375a)2 

- o )/o  = a 1 ^ — — = . 7778 

^ Cr  Cc  (.0375a)2 


from  Equation  (355) 

A 9.8 


a c/6  3 . T O8^ 


= 1.334  • 10”^  sec/sec^ 


and  from  Equation  (368) 


rs  = 


, 2 

V _ .0375 


.1" 


°a  40- 106 


1.334-10 


-8 


= 7.02  • 10 


-4 


(385) 


(386) 


(387) 


Then  from  Figure  42  (for  Ir=32),  using  the  coordinate  values  from 
Equations  (385)  and  (387)  we  obtain  (by  interpolation  between  the  given 
curves ) 

ATf  % .00084  (388) 
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Figure  42.  Normalized  ranging  loop  timing  jitter  vs.  r 
stationary  terminal,  preliminary  analysis. 


Normalized  ranging  loop  timing  jitter  v 
stationary  terminal,  preliminary  analys 


Then 


, .00084  q „ 

=8.4  msec. 


(389) 


The  frame  length  is  then 


I Tf  = 32  • 8. 4 • 10_J  = .269  sec. 


(390) 


From  Section  IV  (see  Figure  20)  we  note  that,  for  stability,  the 
frame  length  must  exceed  the  round  trip  path  delay  which,  for  a 
synchronous  satellite,  varies  from  about  .239  to  .279  sec.  as  the 
terminal  moves  from  the  sub-satellite  point  to  the  horizon.  Conse- 
quently the  above  value  might  be  considered  marginal. 

We  note,  from  Figure  43,  that  if  we  allow  Ir  = 64,  rather  than 
the  original  value  of  32  as  given  above,  the  result  is 


ATf  = .00051 
Tf  = 5.1  msec 


I T^  = .326  sec 


(391) 

(392) 


(393) 


which  is  adequate  for  all  cases.  Hence  we  choose  these  values. 

The  received  ranging  bit  energy  to  noise  density  ratio  may  be 
found,  if  desired,  from  Equation  (359)  and  the  given  value  = .0375  , 
as 


E/Nq  = 23.37  dB 


(394) 


and  if  we  choose  the  number  of  slots  per  subframe  K=32  then  the  ranging 
pulse  (slot)  length  is 


= Tf/K  = 1.59  • 10'H  sec 


(6375  chips) 


(395) 


and  the  received  power  to  noise  density  ratio  is 


P/N0  = V1  = 61.3  dB  Hz. 


(396) 


The  maximum  clock  loop  error  variance  and  the  Kalman  gain  factors 
may  be  obtained  from  Figures  33  through  36,  as  in  the  previous  section, 
using  the  values  of  *Tf  and  rs  as  obtained  above.  The  results  are 


% O o 

P = 4.6  • 10  c 

(397) 

Ky  = P^2  = 4.4  • 10"2 

(398) 

K|  Tf  = 9.8  • 10~4 

(399) 

Ky  Tf2  = 1.1  • 10‘5 

(400) 

from  which 


9.8  • 1 0-4/ 5 . 1 * 10'3  = .192 
1.1  • 1 0-5/ (5.1  • TO'3)2  = .423 


(sec-1 ) 

(401) 

(sec-2) 

(402) 

and  the  standard  deviation  of  the  clock  loop  timing  error  (maximum)  is 
(P^172  = [4.6-10'2(.0375A)2]1/2  = .0080A  . (403) 

Alternatively,  we  might  have  decided  to  retain  the  value  Ir=32 
in  the  above  example  and  increase  the  permissible  subframe  time  Tf  by 
choosing  Nr=2  (i.e.,  by  averaging  2 ranging  pulses  per  frame).  In  this 
case  by  replacing  o^r2  with  °fr2/Nr  in  Equation  (385)  we  obtain 

Kr2  • V2/Nr)/o?c2  = K278  * (404) 

Using  this  value  as  the  ordinate,  together  with  rs=7.02  • 10~4  (from 
Equation  (387)),  we  find  from  Figure  42 

ATf  * .0011  (405) 


then 


Tf  = .0011/A  = 11  ms 


(406) 
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and  the  frame  length  becomes 


I Tf  = 32-1T10'3  = .352  sec.  (407) 

Now,  if  the  2 ranging  pulses  are  spaced  a subframe  (Tf)  apart  as 
indicated  by  Figure  20  (Section  1V-F),  then  we  must  have 

IrTf  * Tmax  + Tf  (403) 


where 


t + T,  = .279  + .011  = .290  sec  . (409; 

max  f 

Hence  this  alternative  is  also  satisfactory.  The  remaining  quantities 
of  interest  may  then  be  found,  if  desired,  as  in  the  above  example. 

It  can  be  seen,  from  Equation  (229)  and  Figures  42-44,  and  similarly 
for  Equations  (234),  (298)  and  (342)  that  there  is  a definite  limit  on 
the  value  of  measurement  noise  permissible  to  obtain  a given  maximum 
ranging  loop  timing  jitter.  To  obtain  a solution  we  must  have 

Kr  - V2/"r)/otc2  * 0 (4'°> 

from  which 


£r 


< o a r 

er  r 


(411) 


This  definite  limit  did  not  exist  in  the  clock  loop  since  the  Kalman 
filter  used  in  that  loop  essentially  provides  optimum  averaging  over 
all  previous  samples  rather  than  the  discrete  number  of  samples,  Nr, 
used  here.  To  obtain  reasonable  solutions  we  would  expect  to  find  the 
normalized  timing  jitter  to  be  near  unity,  and  almost  certainly  in  the 
range 


lO'2  < (a£r2  - o5r2/Nr)/aCc2  < 102  . (412) 

A similar  solution  could  be  obtained  for  the  stationary  satellite 
case  from  Equation  (234)  and  the  curves  of  Figure  8 using  essentially 
the  method  of  Example  1.  However,  this  solution  will  be  obtained  as 
one  term  in  the  solution  for  the  more  exact  expression,  Equation  (298) 
hence  this  equation  will  now  be  evaluated. 
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C.  Application  of  The  Extended  Up-Link  Timing  Analysis 
1 . Stationary  satellite,  maneuvering  terminals 
If  we  let 

Ar  = Nr  = 1 (413) 

(unity  ranging  loop  gain  constant,  and  a single  ranging  pulse  per 
measurement) 

and 


No  = Nd  = Np  - 1 (414) 

(variance  computed  just  prior  to  receipt  of  another  ranging  pulse, 
and  path  delay  k frame  length) 

where  Np  ^ I = the  number  of  subframes  per  frame 


analogous  to  the  special  case  (essentially  worst  case)  considered  in  the 
preliminary  analysis  for  uplink-timing,  stationary  satellite  case;  then, 
with  appropriate  changes  in  notation  to  facilitate  comparison  with  the 
previous  analysis,  Equation  (298)  reduces  to 


P,(I  -1)  P,(I  ) 
1 r 1 r 


P, (21-1 ) 

r 2 -! V-  (415) 


(where  the  subscripts  r and  c on  the  og's  refer  to  the  ranging  loop  and 
clock  loop,  respectively)  while  the  normalized  form  of  Equation  (234) 
from  the  preliminary  analysis  is 


(416) 


The  first  term  on  the  right  hand  side  of  Equation  (415)  (and  the  only 
term  of  Equation  (416))  is  simply  3 tines  the  normalized  clock  loop 
timing  filtered  error  variance,  or  from  Equations  (121)  and  (127), 

3 times  the  Kalman  gain  factor  Kj. 


1 


The  second  and  third  terms  of  Equation  (415)  may  be  evaluated  using 
Equation  (189)  from  which 


P(M  = [(I-KH)*]4  P 


(417) 


where  (from  Equations  (88),  (91)  and  (120)) 


[ ( I - KH ) $ ] = 


1-K, 


(l-KT)Tf 


0-Kt)  2~ 


■kt  ]-KtTf 


TfKT  "IT 


-KTTf 


p-K: 


t r 


(418) 


the  filtered  error  covariance  matrix  for  the  clock  loop,  P,  is  defined 
as 


P = 


(419) 


and 


p(M  5 


^(0 

P2(») 

p4(‘> 

P3(t) 

p5<«) 

£4(0 

‘p6(M 

(420) 


Note  that  while  P (Equation  (419))  is  symmetric,  P(£)  (Equation  (420)),  in 
general,  is  not.  This  is  indicated  by  the  underscore  notation  in  Equation 
(420)  which  permits  retaining  the  original  subscripts  to  avoid  confusion, 
i.e., 


P.  (*•)  / P_.  (*) , in  general . 


(421) 
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From  Equations  (418),  (419)  and  (420)  there  results 


^(0  = 

all 

h + ai  2 

P2  + 

where  the  a^  . are 

el  ements 

of 

an 

a12 

al  3 

[a]  5 

a21 

a22 

a23 

J31 

a32 

a33 

*13  P4 


C(I-kH)*r 


(422) 


(423) 


Unfortunately,  the  above  equations  are  not  convenient  for  calculation 
since  the  matrix  [(I-kH)<f]  is  an  explicit  function  of  Tf  (see  Equation 
(418)).  It  can  be  shown  that  the  normalized  error  variance  (Equation 
(415))  is,  in  fact,  independent  of  Tf;  hence  we  desire  an  alternate 
formulation  in  which  this  pseudo  dependence  upon  Tf  does  not  appear.  One 
way  to  do  this  is  to  redefine  the  system  state  variables  and  the  resultant 
functions  in  dimensionless  (normalized)  form;  i.e..  Equations  ( 83 ) - (91 ) , 
etc.  are  replaced  by 


zVn 


and 


(424) 


where 


Z 


(n) 


Vf 

V? 


0 1 1 

0 0 P 


(425) 


(426) 


(427) 
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K'  = 


KtTf 

L_KTTf 


H‘  = H = [1  0 0] 


r1  = r = 


0 

0 

1 


and 


’(n)  C(n) 


U(n)  U(n)  Tf 


(428) 


(429) 


(430) 


(431) 


(432) 


in  terms  of  this  new  set  of  dimensionless  variables  we  can  obtain 

MO  p. 


_ . f’l  P2Tf 

o2  " bll  ^ + b12  -y.  + b. 


P 

Tf 


T 

where  the  b..  are  elements  of  the  dimensionless  matrix 


(433) 


[b]  h 


1 1 


°21 
L> 


b12 

b13 

b22 

b23 

b32 

b33 

= [ ( I-K' H)4 1 ]* 


(434) 


and 
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[(I-K'H)*']  = 


1-Kt  1-Kt 


1-K, 


-K^-Tf  1 - KjTf  1- 


KfTf 


T2 

2 2 f 

-ICpf  -K”T^  p-K^  ~ 


Then  from  Equations  (120)  and  (433) 

Pi  («•) 

~~2 ” = bll  h + b12  KTTf  + b13  KTTf 


/..t2 


(436) 


The  normalized  gains  Kj,  KjTf,  and  KjT(p  appearing  in  Equations  (435)  and 
(436)  are  directly  obtained  from  the  solution  of  the  clock  loop  filter, 
Equations  (120)  to  (129). 

Curves  of  the  normalized  transmit  timing  error  variance  for  the 
stationary  satellite-moving  terminal  case  versus  the  parameter  rs  as 
obtained  from  Equation  (415)  are  given  in  Figures  45-48  for  different 
values  of  Ir,  and  a comparison  of  tin's  solution  v/ith  the  approximate 
solution  (Equation  (416))  is  given  in  Figure  49.  It  is  seen  that  the 
approximation  is  good  for  higher  values  of  ATf  (lower  maneuver  cor- 
relation) and  for  lower  values  of  rs  (lower  measurement  noise  to  maneuver 
noise  ratio)  with  some  deterioration  elsewhere.  Over  the  range  of 
parameters  investigated  (as  shown  by  the  curves)  the  contribution  from 
the  second  term  of  Equation  (415)  was  found  to  be  at  least  an  order  of 
magnitude  less  than  that  from  the  1st  term  and  hence  could  generally  be 
neglected  with  little  error.  The  third  term,  hov/ever,  becomes  significant 
for  higher  values  of  r^  and/or  lower  values  of  ATf.  Both  the  second  and 
third  terms  alternate  in  sign  thus  producing  the  "ripples"  evident  in  the 
curves  of  Figures  46-49. 

Example  5 

Given  the  system  parameters  of  Example  4 we  have  from  Equation  (385) 


2 2 

a - a, 

er  f’r 


^c 


= .7778 


(437) 


= 64 


and  from  Equation  (387) 


rs  = 7.02  • TO'4  (438) 

Then  from  Figure  46  for  the  stationary  satellite  case  (with  Ir  = 32) 

ATf  = .00500  (439) 


and  using  the  parameters  A = .1  (sec~^)  and  I = 32: 

T^r  = 50.0  ms  (subframe  time)  (440) 

and 

Tp  = 1^1^  = 1.60  sec  (frame  time).  (441) 

A similar  solution  using  curves  generated  from  the  approximate  expression 
(Equation  (416))  yields  essentially  the  same  results  as  above  (Tf=49.9  ms). 

The  effects  of  changes  in  the  various  parameters  upon  the  solution 
for  the  maximum  subframe  time  Tf  (to  insure  a timing  jitter  o£r  <_  the 
specified  value)  are  shown  by  the  families  of  curves  in  Figures  50  to  55. 
These  were  obtained  by  computer  iteration  converging  to  specified  values 
of  normalized  ranging  loop  timing  jitter  in  a manner  analogous  to  that 
used  previously  in  the  study  of  clock  loop  timing  jitter.  As  noted  in 
the  above  examples,  the  minimum  usable  subframe  time  is  determined  by  the 
requirement  that  the  frame  length  Tp  = I rT f must  be  greater  than  the 
round  trip  path  delay.  The  limiting  value  Tp  = 0.3  sec  (applicable  to  a 
synchronous  satellite  transponder  - with  some  safety  factor)  is  also  shown 
in  these  figures.  Usable  values  of  Tf/A  for  system  design  then  lie  below 
the  applicable  curve  of  each  family,  and  above  this  minimum  value. 

As  an  example,  for  a system  with  parameters  as  given  in  Figure  50, 
with  a <_  ,05a  (for  no^signi ficant  increase  in  error  rate)  we  could 
choose  Tf/A  = 3.75  • 10b  (the  minimum  value)  and  thus  accommodate  terminal 
accelerations  considerably  greater  than  lOg.  However,  if  we  desire  to 
limit  the  maximum  error  standard  deviation  aer  to  .0425a  in  order  to  obtain 
an  improved  error  rate  (other  parameters  remaining  unchanged),  then  the 
maximum  permissible  terminal  acceleration  will  be  about  2g.  Note  that 
this  lower  limit  (IrTf  = 0.3  sec)  does  not  apply  to  airborne  relay  systems 
or  other  similar  applications  where  the  propagation  delay  is  short.  In 
such  cases,  much  shorter  frame  lengths  could  be  used,  if  desired,  to 
permit  higher  values  of  terminal  acceleration  and/or  lower  error  rates,  or 
to  compensate  for  desired  changes  in  other  system  parameters.  As  is 
evident  from  these  figures  the  path  lengtii  delay  limitation  (Tf  .3  sec) 
is  generally  not  of  great  significance  over  the  range  of  values  studied 
for  case  of  the  stationary  satellite  relay. 
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2.  Stationary  terminal,  maneuvering  satellite 

We  will  again  consider  essentially  the  "worst  case"  condition 
analogous  to  that  used  in  the  previous  analyses  by  choosing 

N = N.  = Np  , (442) 


where  Np  = I . 

Also  let 

aa  Tf  + To  <<:  Tf  • (443) 

Then,  with  appropriate  change  of  notation  as  used  in  the  previous  section, 
the  normalized  ranging  loop  timing  error  variance  (from  Equation  (342)) 
can  be  written  as 


2 Rj,  * (N  ) 
^lp  £rn  0 0 


(444) 


where  the  subscript  t denotes  the  stationary  terminal  case. 

The  1st  term  on  the  right  hand  side  of  Equation  (444)  is  simply  the 
normalized  timing  error  variance  previously  obtained  for  the  case  of  the 
stationary  satellite  (Equation  (415)),  while  the  remaining  terms  are 
obtained  from  Equations  (336)  and  (345),  respectively  (utilizing 
Equations  (442)  and  (443)  as: 


°o  (V  ? 

q1p  ° = 4(1 -p; 


V2 

'<2  Ir-D  Sp2(Ir-l)  - 2 k|0  S0(k)[Sp(Ir-l>-Sp(k>] 


Rj(0)Tj 


2Ir-2 


+ Cf  <(2I  -1)  — — + 2 l (21-1 -*)  E 

a 1=1 


Rj(OTj 


T 


C I 


y 


'1-2  R=f„(i)T 


T u ' f 


[S(l-l)(2I-l-l)-F(i)]>  (445) 


P r 
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and 


-2  R'v  -v  (N  ) 
q.  e o o 

^ 1 n rn 


21  -2  n(4)T 


^ Sp(Ir-l)  I 


£ u'  ' ’f 
cn 


11=1-1  o. 


where 


I -2 
r 


S U) 


4=0 


*2  u<‘)i  u<irw)T? 


cn 


cn 


+ 2 C 


21-2  j 

r e i 
v cn  e 

L 2 

4=1  -1  c/ 

r £ 


(*)T 


>r  (t)Tf 

v e cn 

L 2 

4=1  a/ 


2Ir-l  Rt  -v 

V 
L 

4=1  + 1 
r 


T „ (On 

e£cn 


(446) 


C = 


Ir-1 

( 1 +p  ) ( 1 -p  r ) 
2 ( 1 -P  ) 


(447) 


4a  ; 


r - 


°aTf 


(448) 


and  Sp(k)  and  Fs(4)  are  as  previously  defined  by  Equations  (316)  and 
(339),  respectively. 
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To  evaluate  Equations  (445)  and  (445)  we  need  expressions  for  the 

correlation  functions  Rv  (*),  Rj  „(*).  R%  „(*).  R*  V f*) , and 
* % / \ g ' e « ccnu  Grn  >e 

RT..ecn  ' Froni  the  definin9  equations  (Equations  (337),  (338),  (346), 

(347),  and  (348))  and  Equations  ( 1 89 ) - ( 1 91 ) we  obtain: 

Rt  (*•)  = P(,(^)  (445 


Rj  u ( ^ ) = °a-[  W 


u(i)  = a.T(A) 


Rg  T (£)  = M*) 
cn  e 


Rj  £ i^)  = Pq(^) 

e cn  — 


where  oflT(«.)  and  ^(M  were  previously  defined  by  Equation  (192)  and 
the  P..(a)  were  defined  by  Equation  (420). 

Finally,  using  methods  analogous  to  those  of  the  previous  section 
there  results 


Ry  (®-)Ty  n y2  P J -r4 

e . h Vf  P5Tf  + . P6Tf 

"2  b31  7 b32  7 b33  7 

5 °C  °£ 

c ^C  ’C  cc 


Rf  u^Tf  -33 ^r 


b,,(£)  (1-p2)  , t > 1 


0 , n < 0 


t 


FK  U(*)T‘ 
cn 


■ 


k,3iF> 


£ > 1 


(456) 


and 


cn  e 


Rt  % (on 

e cn  f 


= b 


11 


i < 0 


- 2 - 3 

P4Tf  , P5Tf 

2 + b12  2 + b13 


P 

V f 
1 


(457) 


P1  P2Tf  Vf 

b31  2 + b32  2 + b33  2 

°C 

c c ^c 


(458) 


where  the  b-jj  are  as  defined  in  Equation  (434),  the  b_-j  ,•  are  defined  as  the 
corresponding  elements  of  the  dimensionless  matrix 


[b]  = 


-11 

-12 

-13 

^21 

^22 

^2  3 

-31 

-32 

-33 

= [(I-K'M)*']*"1 


(459) 


The  normalized  clock  loop  filtered  error  covariance  terms  are  obtained 
from  Equations  ( 1 20 ) - ( 1 29 ) , ( 160)- ( 165) , and  (448)  as 


„ _ JL 
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(462) 
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(464) 


(465) 


where  A,  B,  and  D are  given  by  Equations 
(101),  (100)  and  (97), 


F 


E 

and 

C 


F+l  p BD 

- fTpRTVAT 

Zf  (jf^  - 0)  + 2"  (3E-F-1  ) 


( 1 24 ) - ( 1 26 ) and,  via  Equations 


(466) 


(467) 

(468) 


Curves  of  the  normalized  transmit  timing  error  variance  for  the 
stationary-terminal  maneuveri nq-satel 1 i te  case  versus  the  parameter  r5 
as  obtained  from  Equation  (444)  are  given  in  Figures  56  and  57  for  values 
of  Ir  equal  to  32  and  64  subframes/frame,  respectively.  The  1st  term  of 
the  variance  expression  (Equation  (444))  is  dominant  throughout  most  of 
the  range  of  the  curves.  The  second  term  is  dominant,  however,  for  small 
values  of  rc,  but  decreases  rapidly  with  increasing  rs  thus  producing 
the  "breaks^'  in  slope  of  the  curves  which  are  evident  for  larger  values 
of  ATf.  The  third  term  is  of  alternating  sign  thus  contributing  to  the 
ripples  evident  in  these  curves  for  lower  values  of  rs.  This  term  also 
decreases  with  increasing  rs  becoming  generally  negligible  for  rs  .01. 
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Figure  57.  Normalized  ranging  loop  timing  jitter 


The  ripple  for  higher  values  of  rs  is  present  in  the  first  term  (see 
Equation  (415)  and  the  associated  discussion  of  sub-section  1.  above). 


A comparison  between  this  solution  for  the  ranging  loop  timing  error 
variance  (stationary  terminal)  and  the  approximate  solution  (from 
Equation  (229))  is  given  in  Figure  58.  On  the  average,  the  approximate 
solution  agrees  well  with  the  more  exact  solution  but  the  differences 
may  be  significant  for  typical  values  of  system  parameters  (see  Example  6 
below).  The  approximation  should  be  adequate  for  preliminary  system 
design,  however,  and  offers  a considerable  reduction  in  computation  time, 
particularly  for  large  values  of  Ir. 

Example  6 

Using  the  system  parameters  of  the  previous  example 


2 2 

a -a 

Ert  ar 


= .7778 


and 


rs  = 7,02 


10 


-4 


from  Figure  56  (for  Ir=32)  we  obtain 


XTf  = .00124  . 


(469) 


(470) 


(471) 


Then 


.00124 

— rr 


12.4  ms 


(472) 


and 


IrTf  = 0.397  sec  (473) 

System  parameters  for  this  example  are  compared  to  those  obtained 
for  Examples  4 and  5 in  Table  3. 

Effects  of  changes  in  parameters  upon  the  maximum  usable  subframe 
time  for  the  case  of  a stationary- terminal  maneuvering- satellite  relay 
are  shown  by  the  families  of  curves.  Figures  59  to  64.  These  curves  were 
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Figure  58.  Comparison  of  extended  and  approximate  solutions  for  the 
ranging  loop  timing  jitter,  stationary  terminal. 


Table  3 

SYSTEM  PARAMETERS  FOR  RANGING  LOOP  SYNCHRONIZATION 
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Figure  61.  Subframe  length  vs.  maximum  acceleration  for  different 
maneuver  probabilities,  stationary  terminal. 


obtained  via  an  iteration  technique  using  Equation  (444).  As  for  the 
analogous  curves  (Figures  50  to  55)  of  the  previous  section  (stationary- 
satellite  maneuvering-terminal ) the  maximum  usable  subframe  time  is  given 
by  the  appropriate  curve  of  each  family  while  the  lower  limit  is  set  by 
the  round  trip  path  delay  (t  0.3  sec,  for  a synchronous  satellite). 


It  is  evident  from  these  figures  that  the  maximum  usable  subframe 
time  for  the  maneuvering  satellite  case  is  considerably  less  than  that 
for  the  corresponding  case  of  the  stationary  satellite  (as  would  be 
expected).  Hence,  since  the  lower  bound  (Tp  % .3  sec)  remains  fixed,  the 
range  of  parameters  within  the  permissible  operating  range  is  corres- 
pondingly reduced.  For  example,  from  Figure  59,  if  the  minimum  permissible 
value  for  the  subframe  time  (Tf/A  % 3. 75* 10^)  is  chosen,  then  a maximum 
satellite  relay  acceleration  of  approximately  2g.  can  be  accommodated  in 
order  to  maintain  a < ,05a,  as  compared  to  considerably  more  than  lOg. 
ert  — 

terminal  acceleration  for  the  corresponding  case  with  a stationary 
satell ite. 


As  noted  in  the  previous  section  the  lower  bound  on  subframe  time  is 
not  applicable  for  essentially  earth-bound  relays  where  propagation  delay 
becomes  negligible.  Also,  as  noted  previously,  the  degradation  in  per- 
formance for  the  case  of  the  maneuvering  satellite  relative  to  that  of 
the  stationary  satellite  arises  because  changes  in  satellite  motion  are 
not  known  at  the  terminal  (where  timing  measurements  and  corrections  are 
made)  until  after  the  (one-way)  propagation  delay.  Hence,  for  an  airborne 
relay  where  propagation  delay  is  negligible,  applicable  system  parameters 
must  approach  those  of  the  previous  section  (i.e.,  the  stationary 
satellite  case)  regardless  of  whether  the  terminal  or  the  relay  (or  both) 
are  maneuvering.  Thus,  for  an  atmospheric  airborne  relay,  considerably 
greater  maneuverability  is  permissible  than  for  a satellite  relay  while 
maintaining  the  same  transmit  timing  error  variance.  This  is  indeed 
fortunate  since  in  the  forseeable  future  airborne  relays  appear  much 
more  likely  to  be  highly  maneuverable  than  do  satellites. 
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SECTION  VII 
SIMULATION  STUDIES 


This  section  presents  typical  examples  of  results  obtained  from 
the  simulation  studies  performed  to  investigate  the  clock  loop  transient 
response  and  the  filter  performance  under  conditions  of  no  terminal 
motion  (zero  X-vector)  to  check  the  validity  of  the  Kalman  filter  model. 

A.  Transient  Response 

The  transient  response  of  the  clock  loop  filter  derived  in  Example 
1 of  Section  IV  is  given  in  Figures  65  and  66.  The  parameters  for  this 
filter  were  given  previously  in  Table  2.  At  time  t.=0  (0  iterations), 
when  the  filter  is  first  turned  on,  there  is  a 1/2  chip  delay  in  the 
system  clock  waveform.  The  filter,  starting  from  zero  delay,  proceeds  to 
acquire  the  delay  as  shown  in  Figure  65a.  Measurement  noise  (af=. 05a)  is 
included  in  the  response  shown  in  this  figure,  hence  this  curve ’represents 
the  actual  transient  performance  of  the  filter,  at  each  iteration  (i.e., 
following  the  processing  of  each  data  sample),  for  parameters  as  specified 
in  table  2.  Figure  65b  shows  the  filter  response  for  the  same  situation 
except  the  actual  measurement  noise  has  been  reduced  to  zero  so  that 
the  filter's  transient  response  is  more  readily  seen.  Figure  66  illustrates 
a similar  situation  except  that  in  this  case  the  Kalman  gains  have 
reached  steady  state  prior  to  the  introduction  of  a 1/2  chip  offset 
in  delay.  (This  is  equivalent  to  replacing  the  Kalman  filter  of  the 
above  example  by  the  corresponding  Wiener  filter.)  The  initial  response 
to  the  offset  is  seen  to  be  somewhat  slower,  but  the  overshoot  is  also 
somewhat  less,  resulting  in  a total  effective  settling  time  not  much 
different  from  that  of  the  previous  example.  Hence,  for  this  system 
there  appears  to  be  little  advantage  in  employing  a Kalman  filter.  The 
simpler  Wiener  filter  would  perform  nearly  as  well. 

An  example  for  which  the  Kalman  filter  with  variable  gains  may 
have  a definite  advantage  over  the  fixed  gain  Wiener  filter  is  shown, 
with  and  without  measurement  noise,  in  Figures  67  and  68.  Parameters 
for  this  filter  (which  approximate  the  parameters  of  the  present  modems 
operating  in  the  lower  rate  format  [6],  but  with  greatly  increased 
measurement  noise)  are  given  in  Table  4.  In  this  case,  the  settling 
time  for  the  Kalman  filter  (Figure  67)  is  considerably  shorter  than  for 
the  Wiener  filter  (Figure  68).  However,  the  overshoot  is  somewhat 
greater  and  the  effects  of  noise  spikes  just  after  turn  on  (when  the 
Kalman  gain  is  high)  are  very  severe,  which  could  cause  problems  with 
lock-up.  Note,  However,  that  the  measurement  noise  variance  used  for 
this  example  is  much  greater  than  would  normally  be  used  in  a practical 
system,  therefore  the  effects  of  noise  spikes  in  a practical  system  would 
be  significantly  smaller.  If  violent  maneuvers  (or  sudden  offsets  in 
the  clock  oscillators)  are  expected  during  normal  operation,  then  some 
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DELAY  (CHIPS)  DELAY  (CHIPS) 

0-0  0.5  1.0  0.0  0.5 


O 


(b) 

Figure  66.  Transient  response  of  the  Kalman  filter  of  Example  1 with 
fixed  steady  state  (Wiener  Filter)  gains. 

(a)  With  noise  (b)  No  measurement  noise 
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Figure  67.  Transient  response  of  the  Kalman  filter  specified  in  Table  4. 
(a)  With  noise  (b)  No  measurement  noise 
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(a) 


o 


(b) 


Figure  68.  Transient  response  of  the  Kalman  filter  specified  in  Table  4 
with  fixed  steady  state  (Wiener  Filter)  gains. 

(a)  With  noise  (b)  No  measurement  noise 
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Table  4 


Filter  Parameters  for  the  Response 

Curves  of  Figures  67  and  68. 

Code  Rate  (1/a) 

140  Kbps 

Max  Timing  Error  Std.  Dev.  (5^)  ^ 

.056A 

Measurement  Noise  Std.  Dev.  (a  ) 

.5a 

Max.  Terminal  Acceleration  (A) 

5 g. 

Maneuver  Rate  (A) 

.2  maneuvers/sec. 

Subframe  Time  (T^) 

.0267  sec. 

E/No 

2.9  dB 

,2 

a 

r.  = — — 

4.80 

°a 

Kj  (steady  state  range  gain) 

.013 

K|  (steady  state  velocity  gain) 

2.97  • 10‘3 

K'j  (steady  state  acceleration  gain) 

2.52  • 10"4 

means  of  re-initializing  the  filter  gains  at  those  times  must  be  incor- 
porated; otherwise,  any  advantages  of  the  Kalman  filter  over  the  steady 
state  (Wiener)  filter  will  be  lost  after  initial  system  turn-on. 

B.  Steady  State  Zero  Motion  Response 


Consider  a Kalman  filter  which  has  been  uniquely  determined  by 
specifying  the  process  and  measurement  models  according  to  Equations 
(4)  and  (5)  of  Section  ill,  i.e.. 


Z , = { Z + r u 
n+1  n n 


X„  = H Z + £ 
n n n 


(474) 

(475) 


? 

The  Kalman  predictor  algorithm  uses  $,  r,  H,  and  the  variances  o and 
to  obtain  a prediction  of  the  state  Z via  the  equation 

2 ..  = * 1 + K [X  - H 2 ] (476) 

n+1  n -n  n n ' ' 


where  Kp  is  the  prediction  gain  vector  (Kp  = ®K(n)  of  Section  m).  At 
the  n+1  stage,  the  error  in  the  prediction  is  given  by 
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(477) 


^n+1  2n+l  " n+1 


Upon  substituting  Equations  (297)  and  (299)  into  (300),  one  obtains 


7 ,,  = i>  z + r u - *z  - K [x  -h2  ] 

n+1  n n n -nL  n nJ 


(478) 


Under  the  constraint  of  zero  motion.  Equations  (300)  and  (301)  reduce 
to 


'n+1 


-n+1 


n+1 


($-K  H)  \ - K £ 
v -n- ' n -n  n 


(479) 

(480) 


The  estimation  error  variance  under  this  constraint  is  given  by 

Vi  - ■ <*-y>  Vs-v>T  * v*2  a,t 


(481) 


since  the  measurement  noise  at  the  n^h  sample  is  uncorrelated  with  the 
estimate  at  the  n^h  sample  (Equations  (475)  and  (476)).  If  a steady- 
state  exists  (as  i£  does  for  the  models  of  interest  here),  then  in  the 
steady-state  Pn+i  = Pn= P and  Kn+l  = Krl=K)  so  Equation  (481)  becomes 


£ = (*-KH)  P (*-KH)T  + Ka?2KT 


(482) 


Equation  (482)  is  linear  in  P and  can  easily  be  solved  by  a computer. 
Figure  69  shows  the  variation  in  the  normalized  standard  deviation  of 

the  prediction  error  as  a function  of  the  parameter  rs  = 

for  various  Mf.  The  solid  lines  are  general  Kalman  filter  results 
(no  constraints)  as  obtained  in  Section  IV,  and  given  here  for  com- 
parison, while  the  dotted  lines  are  the  solutions  under  the  zero  motion 
constraint.  Note  that  the  error  standard  deviation  is  reduced  under  the 
zero  motion  constraint  as  stated  previously,  although  the  filter  is  no 
longer  optimum  for  this  condition. 

The  results  of  Figure  69,  for  zero  motion  constraints,  were 
verified  by  a computer  simulation  in  which  the  filter  is  given  a noise 
input  and  the  error  variance  is  observed.  $,  r,  H,  and  the  prior 
statistics  a ^ , and  aj~  were  specified  and  the  filter  was  iterated  with 
no  input  until  steadyustate  was  achieved.  After  steady  state  was 
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achieved,  the  filter  was  given  an  input  of  Gaussian  noise  with  variance 
a ^ and  iterated  for  an  additional  10,000  iterations.  A running  sum  of 
tHe  squared  prediction  error  was  kept  and  after  the  10,000  iterations 
the  sum  was  divided  by  10,000  to  yield  the  error  variance. 

Figure  70  shows  ^ as  a function  of  r^  for  various  Mf.  For 

reference  purposes  the  solid  lines  show  the  calculated  results  while  the 
simulation  results  are  shown  by  the  plotted  points.  In  general,  there 
is  close  agreement  between  the  theoretical  and  the  simulated  results. 

The  differences  occur  for  large  values  of  rs  and  for  small  values  of 
ATf.  This  combination  of  parameters  yields  a filter  with  relatively 
slow  response.  Hence,  simulations  significantly  longer  than  10,000 
iterations  may  be  required  to  achieve  closer  agreement  with  theory  in 
these  cases. 
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SECTION  VIII 
CONCLUSIONS 


Optimum  linear  recursive  (Kalman)  filter  techniques  have  been 
applied  to  the  problem  of  system  synchronization  for  a bit  synchronous 
TDMA  satellite  or  airborne  relay  communications  system  employing  two 
sampled  data  delay  lock  loops  at  each  terminal  to  establish  receive  clock 
and  transmit  timing  relative  to  a system  clock  at  the  satellite. 

A system  containing  a Kalman  filter  in  the  receive  clock  loop  with 
open  loop  coupling  to  a ranging  loop  containing  a simple  averaging  filter 
was  analyzed.  It  was  shown  that  sufficient  timing  accuracy  should  be 
achievable  to  permit  reliable  operation  at  a code  chip  rate  of  40  Mbit/s 
or  greater  for  systems  containing  highly  maneuverable  terminals  and/or 
relays.  For  the  case  of  a maneuvering  satellite  relay  with  stationary 
terminals,  satellite  accelerations  up  to  about  2 g.  could  be  accommodated; 
while  for  a stationary  satellite  and  maneuvering  terminals,  and  also  for 
the  case  of  an  airborne  relay  (where  propagation  delay  effects  on 
determination  of  relay  motion  at  the  terminals  are  negligible),  acceler- 
ations greater  than  10  g.  should  be  possible.  Considerably  higher  code 
rates  and  terminal  maneuverability  (for  given  system  parameters)  could  be 
accommodated  for  down-link  synchronization  only. 

Some  improvements  in  up-link  synchronization,  particularly  with 
regard  to  mean  timing  error,  oscillator  stability  requirements,  and 
lock-up  capability,  may  be  achievable  by  employing  an  optimum  filter 
in  the  ranging  loop  also,  a subject  recommended  for  further  study  if 
higher  code  rates  and  greater  maneuverability  are  desired.  The  tracking 
capability  of  the  ranging  loop,  however,  is  limited  by  the  round  trip 
propagation  delay  and  hence  in  the  case  of  a synchronous  satellite 
relay  a major  portion  of  the  tracking  capability  would  still  need  to  be 
provided  by  open  loop  coupling  from  the  clock  loop  so  that  up-link 
timing  accuracy  must  always  be  somewhat  less  than  the  clock  loop  timing 
accuracy. 

The  analyses  performed  include  development  of  a two-dimensional 
and  a three-dimensional  (augmented)  Kalman  filter  model  for  use  in  the 
clock  loop,  and  for  the  latter  model  a closed  form  steady  state 
solution  was  derived  which  to  our  knowledge  has  not  perviously  been 
obtained.  Steady  state  and  transient  solutions  for  these  models  were 
computed  for  a variety  of  system  parameters  and  are  presented,  together 
with  typical  examples  and  design  curves,  in  a form  particularly  applicable 
to  system  design. 

Simulation  studies  were  also  performed  to  further  investigate 
transient  response,  to  compare  the  performance  of  the  variable  gain 
Kalman  filter  to  that  of  the  fixed  gain  lliener  filter,  and  to  verify 
the  models  and  the  analytical  results  obtained.  Simulation  results  were 
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in  good  agreement  with  analytical  results  in  all  cases,  and  it  was  found 
that  for  high  data  rate  systems  of  this  type,  the  fixed  gain  Wiener  filter 
often  performs  nearly  as  well  as  the  variable  gain  Kalman  filter.  There 
may  be  significant  differences  in  lock-up  capability  however,  and  this 
is  the  subject  of  current  investigation. 

Recommendations  for  future  work  in  this  area  include  design,  bread- 
boarding, and  testing  of  critical  parts  of  the  proposed  TDMA  timing 
system,  particularly  the  clock  loop  Kalman  filter;  followed  by  eventual 
construction  and  testing  of  a complete  timing  system  including  both  the 
clock  loop  and  ranging  loop.  (An  incremental  time  base  shifter  appli- 
cable to  this  system  for  local  clock  timing  corrections  is  currently 
being  developed  under  this  contract.)  Current  theoretical  studies 
involving  lock-up  characteristics  and  transient  response  should  be 
continued  and  an  investigation  of  the  possible  advantages  (primarily  in 
lock-up  and  transient  response)  of  incorporating  an  optimum  filter  in 
the  ranging  loop  is  also  suggested. 
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